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Preface 


The purpose of this publication is to describe the 
FORTRAN library subprograms and their use in either 
a FORTRAN or an assembler language program. The 
body of the publication describes the mathematical 
subprograms (which perform computations) and the 
service subprograms (which perform testing and 
utility functions). This information is intended pri- 
marily for the FORTRAN programmer; Appendix E is 
intended for the assembler language programmer. 
Additional appendixes contain algorithms (the method 
by which a mathematical function is computed ), per- 
formance statistics, descriptions of interruption and 
error procedures, storage estimates, and sample stor- 
age printouts. 

The reader should be familiar with one of the follow- 
ing publications: 


IBM System/360 FORTRAN IV Language, Form 
C28-6515 


IBM System/360 Basic FORTRAN IV Language, 
Form C28-6629 


IBM System/360 Operating System: Assembler Lan- 
guage, Form C28-6514 


Firtu Epirion (October 1968 ) 


IBM System/360 Model 44 Programming System: 
Assembler Language, Form C28-6811 


IBM System/360 Disk and Tape Operating Systems: 
Assembler Language, Form C24-3414 


In addition, references are made within this publica- 
tion to information contained in the following publi- 
cations: 


IBM System/360 Principles of Operation, Form 
A22-6521 


IBM System/360 Operating Systems Supervisor and 
Data Management Macro-Instructions, Form C28-6647 


IBM System/360 Model 44 Programming System: 
Guide to System Use, Form C28-6812 


Standard mathematical notation is used in this publi- 
cation. The reader is expected to be familiar with this 
notation and with common mathematical terminology. 

Appendix E provides information about the use of 
these library routines by assembler language program- 
mers. Corresponding information for pos is in the 
publication IBM System/360 Disk Operating System, 


_ FORTRAN IV Programmer's Guide, Form C28-6397. 


This is a revision of, and makes obsolete, C28-6596-3. Additions have been made to support 
the IBM System/360 Disk Operating System. All changes to the text or tables are indicated by 
a vertical line to the left of the change; revised illustrations are denoted by a ®@ to the left 


of the caption. 


Specifications contained herein are periodically subject to change; any such changes will be 
reported in subsequent revisions or Technical Newsletters. Before using this publication in 
connection with the operation of IBM systems, refer to the latest SRL Newsletter, Form 


N20-0360 for the editions that are applicable and current. 


A form for readers’ comments is provided at the back of this publication. If the form has 
been removed, comments may be addressed to IBM Corporation, Programming Systems 
Publications, 1271 Avenue of the Americas, New York, N.Y. 10020. 
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The ForTRAN Iv library for the System/360 Operating 
System, the Model 44 Programming System and the 
Disk Operating System comprises two types of relo- 
catable subprograms: mathematical subprograms and 
service subprograms. The mathematical subprograms 
correspond to a subprogram defined by a FUNCTION 
statement in a FORTRAN source module. These sub- 
programs always return one answer (function value) 
to the calling module. The service subprograms cor- 
respond to a subprogram defined by a SUBROUTINE 
statement in a FORTRAN source module. These sub- 
programs may or may not return a value to the calling 
module. 

Calls to the library subprograms are either at the 
programmer's request or in response to program re- 
quirements. Under the System/360 Operating System, 
all calls are processed by the linkage editor, which 
takes the subprograms from the library. The library 
subprograms are then combined by the linkage editor 
with the calling module (either an object or a load 
module) into another load module which is ready for 
execution. Under the Model 44 Programming System 
| and the Disk Operating System, the linkage editor 
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takes the subprograms from the library and combines 
them with the calling module into an executable phase. 


The library subprograms may be called in either a 
FORTRAN or an assembler language program. The next 
two sections of this publication contain calling infor- 
mation for the FORTRAN programmer; Appendix E con- 
tains calling information for the assembler language 
programmer. 


Subprogram Names 


Except for a three-character prefix, the names of sub- 
programs described in this publication are the same 
for all systems. The prefix is uniformly mic for the 
Operating System library, Boa for the Model 44 Pro- 
gramming System, and wr for the Disk Operating 
System. These characters occupy the portion of the 
subprogram name shown as xxx in the text and tables. 
Thus, the subprogram named here as xxxcstoc has 
the name IHccsLoc under the Operating System, 
BOACSLOG under the Model 44 Programming System, 
and iLFcsLoc under the Disk Operating System. 
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Mathematical Subprograms 


The mathematical subprograms supplied in the 
FORTRAN library perform computations frequently 
needed by the applications programmer. The mathe- 
matical subprograms are called in two ways: explicitly, 
when the programmer includes the appropriate entry 
name in a source language statement (see Table 1); 
and implicitly, when certain notation (e.g., raising a 
number to a power) appears within a source language 
statement (see Table 6). 

The following text describes the individual mathe- 
matical subprograms and explains their use in a 
FORTRAN program. Detailed information about the 
actual method of computation used in each subpro- 
gram, the performance of the subprogram, interrup- 
tion and error procedures, and storage estimates can 
be found in the appendixes of this publication. 


Explicitly Called Subprograms 


Each explicitly called subprogram performs one or 
more mathematical functions. Each mathematical 
function is identified by a unique entry name that 
differs from the name of the subprogram. 

A subprogram is called whenever the appropriate 
entry name is included in a rorTRAN arithmetic expres- 
_sion. The programmer must also supply one or more 
arguments. These arguments follow the entry name 
and are separated by commas; the list of arguments is 
enclosed in parentheses. : 

For example, the source statement 


RESULT = SIN (RADIAN) 


causes the rHcsscn subprogram to be called. The sine 
of the value in RADIAN is computed and the function 
value is stored in RESULT. 

In the following example, the nicssgrt subprogram 
is called to compute the square root of the value in 
AMNT. The function value is then added to the value 
in stock and the result is stored in ANs. 


ANS = STOCK + SQRT (AMNT) 


The explicitly called subprograms are described in 
the tables that make up the rest of this section. 
These tables show the general function, subprogram 
name, the ForTRAN library that contains the subpro- 
gram, definition, entry name(s), argument information, 
type of function value returned, and assembler re- 
quirements. The following column headings are used 
in the tables: 

General Function: This column states the nature of 
the computation performed by the subprogram. 
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Subprogram Name: This column gives the module 
name of the subprogram. The first three characters 
are uniformly mc for subprograms in the System/360 
Operating System library, Boa for those in the Model 
44 Programming System library, and mr for those in 
the Disk Operating System library. These three char- 
acters occupy the portion of the subprogram name 
that is shown as xxx in the text and tables. 

Subset: This column indicates those subprograms 
that belong to the Basic rorrran rv (£) library. Unless 
otherwise indicated, all subprograms that belong to 
the £ library also belong to the FORTRAN Iv (G, H, MODEL 
44, and pos) library. 

Definition: This column gives a mathematical equa- 
tion that represents the computation. An alternate 
equation is given in those cases where there is another 
way of representing the computation in mathematical 
notation. (For example, the square root can be repre- 
sented either as y= \/x or y = x!/2.) The definition 
for those subprograms that accept complex arguments 
contains the notation z = x,+ Xi. 

Entry Name: This column gives the entry name that 
the programmer must use to call the subprogram. A 
subprogram may have more than one entry name; the 
particular entry name used depends upon the compu- 
tation to be performed. For example, the macsscn sub- 
program has two entry names: cos and sin. If the 
cosine is to be computed, entry name cos is used; if 
the sine is to be computed, entry name sIN is used. 

Argument Number: This column gives the number 
of arguments that the programmer must supply. 

Argument Type: This column describes the mode 
and length of the argument. INTEGER, REAL, and 
COMPLEX represent the type of number; the notation 
*4, *8, and *16 represent the size of the argument in 
storage locations. 

Note: In Basic FORTRAN tv, a real argument corre- 

sponds to the REAL *4 argument, and a double- 

precision argument corresponds to the REAL *8 argu- 

ment. Complex arguments cannot be used with a 

Basic FORTRAN Iv compiler. 

Argument Range: This column gives the valid range 
for an argument. If the argument is not within this 
range, an error message is issued and execution of this 
load module is terminated. Appendix C contains a 
description of the error messages. 

Function Value Type: This column describes the 
type of function value returned by the subprogram. 
The notation used is the same as that used for the 
argument type. 


Table 1. Explicitly Called Mathematical Subprograms 


Specific Function Subprogram Name Entry Name(s) 



































Logarithmic and exponential Common and natural logarithm xxx CLLOG* CDLOG 
subprograms (described in Table 2) xxxCSLOG* CLOG 
xxx_LLOG DLOG, DLOG10 
xxxSLOG ALOG, ALOG10 








xxxCLEXP* 
xxxCSEXP* 
xxx LEXP 

xxxSEXP 


Exponential 





Square root xxxCLSQT CDSQRT 



























xxxCSSQT* CSQRT 

xxxLSQRT DSQRT 

xxxSSQRT SQRT 
Trigonometric subprograms Arcsine and arccosine xxx LASCN* DARSIN, DARCOS 
(described in Table 3) xxxSASCN* ARSIN, ARCOS 

Arctangent THCLATAN DATAN 
xxx LATN2* DATAN, DATAN2 
IHCSATAN ATAN 


xxxSATN2* ATAN, ATAN2 


CDSIN, CDCOS 

















xxxCLSCN* 


Sine and cosine 


xxxCSSCN* CSIN, CCOS 
xxxLSCN DSIN, DCOS 
xxxSSCN SIN, COS 









Tangent and cotangent xxx LTNCT* DTAN, DCOTAN 
xxxSTNCT* TAN, COTAN 

Hyperbolic sine and cosine xxx LSCNH* DSINH, DCOSH 
xxxSSCNH* 


SINH, COSH 
Hyperbolic tangent xxx LTANH DTANH 
xxxSTANH 


TANH 
Absolute value xxxCLABS* 
xxxCSABS* 


CDABS 
CABS ; 
xxxLERF* 
xxxSERF* 






Hyperbolic function subprograms 
( described in Table 4) 






Miscellaneous subprograms 
(described in Table 5) 






Error function 







DERF, DERFC 
ERF, ERFC 


DGAMMA, DLGAMA 
GAMMA, ALGAMA 








Gamma and log-gamma xxxLGAMA* 


xxxSGAMA* 
















Maximum and minimum value xxx FMAXD DMAX]1, DMIN1 
xxxF MAXI AMAXO, AMINO, MAX0, MINO 
xxxFMAXR AMAXI1, AMIN1, MAX], MIN1 







MOD 
AMOD, DMOD 


AINT 
IDINT, INT 


Modular arithmetic IHCFMODI 
IHCFMODR 

Truncation IHCFAINT 
IHCFIFIX 





*Not available in FORTRAN IV (E) 


Assembler Requirements: This column gives the reg- 2, Floating-point register 4 is used for intermediate 
isters used by the subprogram and the minimum save computation. 
area that the assembler language programmer must 3. The save area must be at least nine full-words in 
supply. For example, the assembler requirements for length. 
the eerie oar ae Detailed information for the assembler language pro- 
save area OF grammer is given in Appendix E. 
This information specifies that: Note: In the following tables, the approximate value 
1. The function value is found in floating-point reg- _ of 218 » 7 is .82354966406249996D + 06; the approximate 
isters 0 and 2. value of 2°° ¢ x is .35371188737802239D + 16. 
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e Table 2. Logarithmic and Exponential Subprograms 





rn t(s) Function | Assembler 
General Subprogram | Sub- aaah e Value Require- 
Function Name set Definition Type’ ments 


| Common xxxCLLOG No y=Pv loge (z) CDLOG | 1 | complex | z40 + 0i complex | registers 
and natural See Note 2 *16 See Note 3 *16 0,2 
logarithm save area 
8F 
| xxx CSLOG N =pv log. (z) CLOG complex | z+#0 + 0i complex | registers 
See Note 2 *8 See Note 3 *8 0,2 
save area 
8F 
xxxSLOG 
N 
N 
Y 
Y 


ALOG 
ALOG10 
‘ CDEXP 
See Note 4 
4 CEXP 
See Note 4 


>0 
>0 
real *4 >0 | real *4 registers 
. 0 (2) 
save area 
5Ft 
real *4 x > 0 real *4 registers 
0 (2) 
save area 
5Ft 
complex | x: S 174.673 complex | registers 
*16 [xo] < (2% +m) | *1 0, 2 
save area 
8F 
complex | x: S 174.673 complex | registers 
*8 |xo] << (2°*m) | *8 0, 2 
save area 
8F 
DEXP real *8 x S 174.673 real *8 registers 
0 (2) 
save area 
OF 
y=e EXP 1 real *4 x S 174.673 rea] *4 register 0 
save area 
12F 
CDSQRT}| 1 complex | any complex complex | registers 
*16 argument *16 0,2 (4) 
See Note 3 save area 
OF 
CSQRT 1 | complex | any complex complex | registers 
*8 argument *8 0,2 (4) 
See Note 3 save area 
OF 
xxxLSQRT x DSQRT 1 | real*8 | x20 real *8 | registers 
= 0 (2, 4) 
save area 
5Ft 
xxxSSQRT Yes 1 | real *4 x=0 real *4 registers 
0 (4) 
save area 
5Ft 
NOTEs: 


. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 
the REAL *8 argument. = 
. PV = principal value. The answer given is from that point where the imaginary part (yz) lies between —a and +7. More specifi- 
cally: —m < ye S 7, unless x1 <0 and x2 = —O, in which case, ye = —7. 
. Floating-point overflow can occur. 
. Where z is a complex number of the form x: + xzi. 
+For OS/360, seven words are needed. 


DLOGI10 real *8 real *8 registers 
0 (2) 
save area 
OF 
xxxCSEXP 


_ 


xxxLLOG y=logex or DLOG real *8 x real *8 registers 
0 (2) 
save area 
OF 
x 
Exponential | xxxCLEXP = 


XxXxLEXP a 
xxxSEXP a 


Square root | xxxCLSQT 


xxxCSSQT 


Oo 
es 
es 

Oo 
Oo 
es 
es 

No 

No 

Yes 





Table 3. Trigonometric Subprograms 























Assembler 
Require- 
ments 


Function 
Value 
Type’ 









General Subprogram 


Function 


Entry Argument(s) 
Definition Name No. Type’ 


Range 
y = arcsin (x) DARSIN 
y = arccos (x) wee 
No y = arcsin (x) ARSIN Ix} 1 
y = arccos (x) ARCOS 
See y = arctan (x) DATAN 1 real *8 
Note 
2 
See y = arctan (x) DATAN 1 real *8 
Note 
2 
" _ 
y = arctan a 































real 
*8 (in 
radians ) 


registers 
0 (2,4) 
save area 
13F 


xxx LASCN 





Arcsine and. 
arccosine 







real 
*8 (in 
radians ) 







registers 
0 (2, 4) 
save area 
13F 







real 
*4 (in 
radians ) 


xxxSASCN registers 
0 (2, 4) 
save area 


10F 



















real 
*4 (in 
radians ) 


registers 
0 (2, 4) 
save area 


10F 






























real 
*8 (in 
radians ) 


any real 
argument 


IHCLATAN registers 
0 (2, 4, 6) 
save area 


5F 


Arctangent 

















any real real registers 
argument *8 (in 0 (2, 4, 6) 
radians) | save area 


5Ft 


xxx LATN2 





















any real real registers 
arguments *8 {in 0 (2, 4, 6) 
(except 0, 0) radians) | save area 










SFt 
















































IHCSATAN | See y = arctan (x) ATAN 1 real *4 any real real. registers 
Note argument *4 (in 0 (2, 4, 6) 
2 radians) | save area 
5F 
xxxSATN2 See y = arctan (x) ATAN 1 real *4 any real real registers 
Note argument *4 (in 0 (2, 4, 6) 
P radians) | save area 
5FFT 

























= arctan *) any real real registers 
Urroeas Xe arguments *4 (in 0 (2, 4, 6) | 
(except 0, 0) radians) | save area 






5Ft 


- 
No CDSIN 1 complex complex 
*16 (in *16 
radians ) 
y = cos (z) CDCOS complex | |x:|<(2°+ 7) complex 
See Note 4 *16 (in |xe| 174.673 | *16 
radians ) 
y = sin (z) CSIN |x1]<C( 2 * or) 
See Note 4 





xxxCLSCN 









y = sin (z) 


See Note 4 


Ixi]<( 950 ba 1) 
|xo| 174.673 


registers 
0, 2 (4) 
save area 
OF 


Sine and 
cosine 





















registers 
0, 2 (4) 
save area 


OF 













xxxCSSCN 












complex 
*8 


1 complex 
*8 (in 
radians ) 


registers 
0, 2 (4) 
Save area 


OF 





































y = cos (z) CCOos 1 complex | |x:|<(2%+ 7) complex | registers 
See Note 4 *8 (in |xe| S$ 174.673 | *8 0, 2 (4) 
radians ) save area 
OF 
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Table 3. Trigonometric Subprograms (Continued ) 


Subprogram | Sub- 
Name set Definition 


xxxLSCN 


xxxSSCN Yes 
- 
xxxLTNCT No y = tan (x) DTAN 1 | real [x|<( 2 + 3) real *8 
*8 (in See Note 3 
radians ) | 
y = cotan (x) DCOTAN| 1 | real |x}<(2 * 3) real *8 
*8 (in See Note 3 
radians ) 
y 


xxxSTNCT 


1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds to 
the REAL *8 argument. . 


2. Instead of the IHCLATAN and IHCSATAN subprograms contained in the FORTRAN IV (E) library, the FORTRAN IV library 
contains the xxxLATN2 and xxxSATN2 subprograms. 


3. The argument for the cotangent functions may not be near a multiple of 7; the argument for the tangent functions may not be near 
an odd multiple of 7/2. 





























Assembler 
Require- 
ments 


Function 
Value 
Type’ 


Argument(s) 









General 
Function 
































Sine and fy Ix] <( 2» ) real *8 AEN 

cosine in ‘ 

(continued ) radians ) save area 
5Ft 















y = cos (x) real [x|<( 2° + 3) real *8 | registers 
*8 (in 0 (2, 4) 
radians ) save area 

5Ft 













real |x| <(2"* + 3) real*4 | registers 

*4 (in 0 (2, 4) 

radians ) save area 
5Ft 










real [x|<( 2% + wr) real*4 | registers 

*4 (in 0 (2,4) 

radians ) save area 
5Ft 












Tangent and registers 

cotangent 0 (2, 4, 6) 
save area 
5FF 











registers 

0 (2, 4, 6) 
Save area 
5Ft 





























TAN 1 | real |x|<(2"* + 3) registers 
*4 (in See Note 3 0 (2, 4) 
radians ) save area 

5F 





















COTAN | 1 | real Ix|<(2'* + 3) registers 
*4 (in See Note 3 0 (2, 4) 
radians ) save area 

5F 


Nores: 


4, Where z is a complex number of the form x +xzi. 
For OS/360, seven words are needed. 
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Table 4. Hyperbolic Function Subprograms 


Sub- 
set 
No 


xxxLSCNH : 
xxxSSCNH 


General 
Function 


Subprogram 


Name Definition 


Hyperbolic 
sine and 
cosine 


Hyperbolic 
tangent 


xxx LTANH 
xxxS TANH 


NOTE: 


Entry 
Name 


- | 


real *4 


Argument(s) 


Range 


|x| <174.673 


|x|<174.673 


|x|<<174.673 


|x| <174.673 


any real 


argument 


any real 
argument 





Function 
Value 
Type’ 


Assembler 
Require- 
ments 


registers 
0 (2,4) 
save area 
OF 


registers 
0 (2,4) 
save area 
OF 


registers 
0 (2,4) 
save area 


OF 


registers 
0 (2, 4) 
save area 
oF 


registers 
0 (2) 
save area 


5F 


registers 
0 (2) 
save area 
5F 


1. In FORTRAN IV (e), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 


to the REAL *8 argument. 


Table 5. Miscellaneous Mathematical Subprograms 


Sub- 

set Definition 
xxxCLABS ‘ y=|z|=(x2+x.?)¥? 
xxxCSABS ‘ y=|z| = (x2+x.?)” 


No 
No 
No 
TT. 
2 ("3 
=—= J e* du 
rs Bi 
y=l-—erf (x) 


General 
Function 


Subprogram 
Name 


Absolute 
value 


Error 
function 





Entry 

Name 
CDABS 
CABS 
DERF 
DERFC 


complex 


* 


1 | complex 
*16 
8 


any complex 
argument 
See Note 2 


any complex 
argument 
See Note 2 


any real 
argument 


any real 
argument 


Function 
Value 





Assembler 
Require- 
ments 


registers 
0, 2 (4) 
save area 


8F 


registers 
0, 2 (4) 
save area 
8F 


registers 
0 (2, 4, 6) 
save area 


11F 


registers 
0 (2, 4, 6) 
save area 


11F 
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Table 5. Miscellaneous Mathematical Subprograms (Continued ) 


Ag t(s) Function } Assembler 
General Subprogram | Sub- Entry SUAS Value Require- 
Function Name set Definition Name No. Type’ ments 


ERF 1 jreal *4 any real real *4 registers 
argument 0 (2, 4, 6) 
save area 
11F 
real *4 any real real *4 registers 
argument 0 (2, 4,6) 
save area 
11F 
real *8 x > 2° and real *8 registers 
x < 57.5744 0 (2, 4, 6) 
save area 
11IF 
real *8 x > 0 and registers 
x <4,2913+ 10% 0 (2, 4, 6) 
save area 
11lF 
real *4 x > 2°? and real *4 registers 
x < 57.5744 0 (2, 4, 6) 
save area 
11IF 
real *4 x > Oand real *4 registers 
x <4,2913 +10" 0 (2, 4, 6) 
save area 
liF 
real *8 any real real *8 register 0 
arguments save area 
OF 
real *8 any real real *8 register 0 
arguments save area 
oF 


Error xxxSERF 
function 
(continued ) 


DGAMMA 


y=loge I’ (x) or DLGAMA 
y=log. f u** e™ du 


vy) 


Gamma and { xxxLGAMA 
log-gamma 
xxxSGAMA 


y=log. I’ (x) or 


foo] 
y=loge f u**e™ du 
0 


y=max (m,... 


Maximum xxx MAXD 
and 

minimum 

values 


IV 
bo 


IIV 
bo 


y=min (x,... 


IV 
bo 


y=max (x,. 


12 


IIV 
bo 


IIV 


iW) 


IV 
bo 


o 
< 


xxx MAXI 


integer | any integer real *4 register 0 
*4 arguments save area 
OF 


*4 
*4 


*4 


any integer 
arguments 


any integer 
arguments 


any integer 
arguments 


¥4 


*4 


register 
See Note 3 
save area 
OF 


register 0 
save area 
OF 


register 
See Note 3 
save area 
oF 





























Assembler 
Require- 
ments 


Function 
Value 
Type’ 


*4 
, 
f 





Table 5. Miscellaneous Mathematical Subprograms (Continued ) 
General 


Subprogram |Sub- Entry Argument(s ) 
Function Name set Definition Name 


Maximum |xxxFMAXR- [Yes |y=max(x,...,xn) | AMAXL =2 |real *4 
and 
Minimum 
Values 
(continued ) MAXI 2 
y=min (x1,...,xXn) | AMINIL 2 
7 | 



















any real 
arguments 


register 0 
save area 
OF 

































IV 


any real 
arguments 


register 
See Note 3 
save area 
OF 





IIV 


any real 
arguments 


register 0 
save area 


oF 

























IV 


any real 
arguments 


register 
See Note 3 
save area 
OF 


Modular IHCFMODI | See | y = x (modulo x2) MOD integer x2 +0 integer register 
arithmetic Note] See Note 5 See Note 6 See Note 3 
4 save area 
OF 





register 0 
save area 
OF 


real *4 X20 real *4 
See Note 6 

real *8 Xe 0 real *8 
See Note 6 

real *4 any real real *4 
argument 








register 0 
save area 
OF 










IHCFMODR | See { y = x: (modulo x) AMOD 
Note} See Note 5 
4 
y = x: (modulo x:) DMOD 
See Note 5 


Truncation | IHCFAINT See AINT 
Note 
4 


IHCFIFIX See IDINT real *8 integer 
Note 
4 


NOTEs: 


1. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument, and a double-precision argument corresponds 
to the REAL *8 argument. 


2. Floating-point overflow can occur. 








y = (signx)*n 
where n is the largest 
integer & |x| 


register 0 
save area 
OF 



















y = (signx)*n 
where n is the largest 
integer S |x| 


register 0 
See Note 3 
save area 

OF 


any real 
argument 




















any real 
argument 


register 
See Note 3 
save area 


OF 








ee 


3. The result is stored in general register 0. 

4, The coding that performs this function is out-of-line in FORTRAN IV (E) and in-line in FORTRAN IV. Out-of-line coding is 
taken from the FORTRAN library by the linkage editor and processed with the calling module. In-line coding is inserted by the 
FORTRAN compiler at the point in the source module where the function is referenced. This means that the in-line functions 
are available in FORTRAN IV by using the appropriate entry name but that they are not part of the library. 


5. The expression x, (modulo xz) is defined as x1 — [=] * xe, where the brackets indicate that an integer is used. The largest in- 
X2 , 


teger whose magnitude does not exceed the magnitude of * is used. The sign of the integer is the same as the sign of = : 
Xe 2 


6. If x: = 0, then the modulus function is mathematically undefined. In addition, a divide exception is recognized and an interrup- 
tion occurs. (A detailed description of the interruption procedure is given in Appendix C.) 
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Implicitly Called Subprograms general function, subprogram name, the FORTRAN 


The implicitly called subprograms perform operations library that contains the subprogram, implicit pone 
required: by the appearance of certain notation in a reference, entry name, argument information, Des? 
pCR RGIS SEATGIIGAE: WWKGAS ALGO ie tO be function value returned, and assembler requirements. 
raised to a power or when multiplication and division The column headed “Implicit eee Refer ene 
of complex numbers are to be performed, the Fortran ——8'V€S 4 sample sour ia ora eae i f he n 
compiler generates the instructions necessary to call a -RORTRAN: SOULE TNO Bue Caney ae a eee 
the appropriate subprogram. For example, if the fol- to be called. The rest of the column headings in Table 
lowing source statement appears in a source module, 6 have the same meaning as those used with the ex- 
enseeeeon plicitly called subprograms. 
aie The action taken within the subprogram depends 
upon the type of base and exponent used. Tables 7 
through 10 show the result of an exponentiation per- 
formed with the different combinations and values of 
The implicitly called subprograms in the ForrrAN __ base and exponent. In these tables, I and J are integers; 
library are described in Table 6. This table shows the _—_A and B are real numbers; C is a complex number. 


where BASE and EXPON are values of the form REAL *4, 
the xxxFRxPR subprogram is called by the FORTRAN 
compiler. 
























Implicit 
Function 







Function 
Value 
Type® 


Table 6. Implicitly Called Mathematical Subprograms 
Subprogram | Sub- Argument(s) 
Name set Reference’ Name 


3 
xxxCLAS No y = 21" Ze CDMPY# 2 } complex *16 
save area 5F 
y = 2/22 CDDVD# 2 | complex *16 registers 0, 2 (4, 6) 
save area 5F 
xxxCSAS No y = 21" Ze CMPY#+ 2 | complex *8 complex *8 registers 0, 2 (4, 6) 
save area 5F 
y = 2/22 CDVD# 2 {complex *8 complex *8 registers 0, 2 (4, 6) 
save area 5F 
xxxFIXPI Yes 54 FIXPI# 2 |i=integer*4| integer *4 register 0 
j = integer *4 See Note 4 
save area 18F 
xxx RXPI Yes y =a**j FRXPI# 2 real *4 real *4 register 0 
j = integer save area 18F 
xxxFDXPI Yes y **; FDXPI# 2 real *8 real *8 register 0 
integer*4 save area 18F 
























General 
Function 


Multiply and 
divide complex 
numbers 


Entry’ Assembler 


Requirements 
registers 0, 2 (4, 6) 










complex *16 





complex *16 









Raise an integer 
to an integral 
power 


gees | 


y 


















Raise a real 
number to an 
integral power 


ll 
i>) 
| 


ll 
=o Jo 


oe | op 
| 
mort tee oy 


a 
a**h FRXPR# 





Raise a real 


xxxF RXPR Yes y= 2 eal *4 real *4 register 0 
number to a eal *4 save area 18F 
real power xxxFDXPD y=a**h | FDXPD# 


eal *8 real *8 register 0 
eal *8 save area 18F 
xxxFCDXI = complex {| complex *16 register 0, 2 
*16 save area 18F 
xxxFCXPI register 0, 2 
.| Notes: 


j = integer 
save area 18F 
1. This is only a representation of a FORTRAN statement; it is not the only way the subprogram may be called. 


2. This name must be used in an assembler language program to call the subprogram; the character # is a part of the name and 
must be included. 


3. In FORTRAN IV (E), a real argument corresponds to the REAL *4 argument and a double precision argument corresponds to 
the REAL *8 argument. 


. The result is stored in general register 0. 












Raise a complex 
number to an 
integral power 


FCDXI# 2 


N 













Zz 





FCXPI# 2 






Z = complex 
*8 


j = integer 
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® Table 7. Exponentiation With Integer Base and Exponent 


ye es <r 
ai I 


Compute the } Function 
function value | value = 1 













I>0 Function value = 1 
ifI=1 

Otherwise, function 

value = 0 



































Function 
value = 0 


Error message 
xxx241T; 
or OA241I 


Compute the | Function 
function value} value = 1 


Error message 
xxx24]T; 
or OA241T 


Function value= — 1 
if = —l and 
if J is an odd 
number 
Function value = 1 
if I = —land 
if J is an even 
number 
Otherwise, function 
value = 0 

























@ Table 8. Exponentiation with Real Base and Integer Exponent 














Exponent (J) 


Compute the {| Function Compute the 
function value} value = 1 function value 


Function Error message Error message 


Base (A) 


A=0 


























value = 0 xxx242I or xxx2492] or 
xxx2431; xxx243I; 
or OA2421 or OA2421 
or OA243I or OA243I 












Compute the 
function value 


Function Compute the 
value = 1 function value 


@ Table 9. Exponentiation with Real Base and Exponent 









Exponent (B) 


ys a 
oO esi kao (eo 


A>O | Compute the | Function Compute the 
function value| value = 1 function value 
A=0 


Function Error message Error message 





















value = 0 xxx244lI or xxx244l or 
xxx245]; xxx245]; 
or OA244I or OA244T 





or OA245I 
A<0 | Error message | Function 
| Xxx253] or value = 1 
xxx263] 


Optional Extended Error Message Facility! 


When the rorrran library has been created by sys- 
tem generation, the extended error message facility 
may be selected. This facility provides the following 
features: 

1, Execution can continue after an error occurs in a 
mathematical function, with a standard result supplied 
as the functional value. 

2. A user-written routine can be supplied to take 
corrective action. 

For OS/360 Full-Language FORTRAN only. 


or OA245] 


Error message 
xxx253I1 or 
xxx263] 







e Table 10. Exponentiation with Complex Base and Integer 


Exponent 


Exponent (J) 


Se BR: 0 
R > 0 and Compute the | Function Compute the 
Ri > 0 function value | value=1 + Oi | function value 
R>Oand | Computethe | Function Compute the 
Ri = 0 function value | value=1 + Oi | function value 
R > O and Compute the [{ Function Compute the 
Ri < 0 function value | value=1 + Oi | function value 
R = Oand Compute the | Function Compute the 
Ri > 0 function value | value=1 + Oi | function value 
R = Oand Function value| Error message | Error message 
Ri = 0 xxx2461 or xxx246I or 
xxx247]1; xxx247]1; 


or OA246I1 or OA2461 
or OA2471 or OA2471 


R = O and Compute the | Function Compute the 

Ri<0 function value | value=] + Oi | function value 

R <0 and Compute the | Function 

Ri >0 function value | value=1 + Oi 
Compute the | Function 
function value | value=1 + 0i 


R < 0 and 

Ri =0 
Compute the | Function 
function value | value=1 + Oi 


Base (C) 


Compute the 
function value 


Compute the 
function value 


Compute the 
function value 


R < 0 and 
Ri<0 





3. Installation-written function subprograms may 
use the FORTRAN extended error message facility for 
the control and printing of error messages comparable 
to those put out by the rorrran library. 

Table 11 describes the function values supplied as 
standard corrective actions and the parameters upon 
which recognition is based when user-supplied correc- 
tive action is chosen. Table 12 shows corrective action 
after interrupt occurrence. 


User-Supplied Corrective Action 


If the user chooses to write his own routine for correc- 
tive action, the user-supplied routine is called as 
follows: 


CALL NAME (IRETCD, IERNO, DATA1, DATAQ,...) 


Value on entry to 
Parameter Description subroutine NAME from NAME 


IRETCD Return code 1 lorO 
IERNO- Error Code Cannot be changed Cannot be changed 


DATAI Data Data in error For IRETCD=1, 

DATA2 upon entry the new values to 
be used in recom- 
puting the function 
are placed in 
DATAI, DATAQ2. 
For IRETCD=0, 
these parameters 
may not be altered. 


Value upon entry 





The subprogram NAME may do the following for 
these parameters: 
e Examine the values in DATA1,DATA2. 
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e If the function is to be recomputed, new data is Action” in Table 11 correspond to pata1,paTa2. When 


placed in DATA1,DATA2. writing the subprogram NAME in FORTRAN, declaration 
e If standard corrective action is desired, rRETCD is must be made for the type of constant being processed. 
changed to 0. For example, for error 246, pata1 (C) is COMPLEX*s8, 


The parameters listed in “User-Supplied Corrective and pataz (1) is INTEGER*4, 


Table 11. Optional Service for Error Situations 



























Options 










Invalid Standard User-Supplied 
FORTRAN Argument Corrective Corrective Action 
Reference Range Action (See Note 1) 






For errors 207-210, see Table 12 
CALL SLITE (I) 












The call is treated as 
a no operation 


I>4 



































































































CALL SLITET (I, J) I>4 J=2 

K=I**J I=0, JS0 K=0 

Y=X**] X=0,IS0 Y=0 

DA=D**I D=0, ISO DA=0 

XA=X**Y X=0, YS0 XA=0 

DA=D**DB D=0, DBS0 DA=0 

CA=C**] C=0+0i, ISO CA=0+0i 

CDA=CD*I C=0+0i, ISO CA=0+0i 

Y=SQRT (X) X<0 Y=|x|*? x 

Y=EXP (X) X> 174.673 Yy=* x 

Y=ALOG (X) X=0 =—* x 
X<0 Y=log |X} Xx 

Y=ALOGIO (X) X=0 =—* xX 
xX<0 Y=logo|X| Xx 

Y=COS (X) |X| 22" Y= V2/2 x 

Y=SIN (X) 

Y=ATAN2 (X, XA) X=0, XA=0 X, XA 

Y=SINH (X) |X|=174.673 Xx 

Y=COSH (X) 

Y=ARSIN (X) |X|>1 Xx 

Y=ARCOS (X) 

Y=TAN (X) |X| = (2%) *x Xx 

Y=COTAN (X) 

Y=TAN (X) X is too close to an Xx 











odd multiple of = 












































Y=COTAN (X) X is too close to a xX 
multiple of x 

DA=DSQRT (D) D<0 DA=|D/” D 

DA=DEXP (D) D>174.673 DA=* D 

DA=DLOG (D) D=0 DA=-* D 
D<0 DA=log]X| D 

DA=DLOGI10 (D) D=0 DA=-—* D 
D<0 DA=logu|X! D 

DA=DSIN (D) [D| 22°" DA= V2/2 D 

DA=DCOS (D) 

DA=DATAN2 (D, DB) D=0, DB=0 D, DB 









Variable Type 









I,J Variables of INTEGER*4 

X, XA, Y Variables of REAL*4 

D, DA, DB Variables of REAL*8 

C,CA Variables of COMPLEX*8 

Z, X1, Xe Complex variables to be given the length of the functioned argument when they appear 
CD Variables of COMPLEX*16 






NOTEs: 

1. The user-supplied answer is obtained by recomputation of the function using the value set by the user routine for the parameters 
listed. (See “User-Supplied Corrective Action.” ) 

2. The largest number that can be represented in floating point is indicated by *. 
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Table 11. Optional Service for Error Situations (Continued ). 





















Invalid Standard User-Supplied 
FORTRAN Argument Corrective Corrective Action 
Reference Range Action (See Note 1) 













266 DA=DSINH (D) |D|=174.673 DA=* D 
DA=DCOSH (D) 
267 DA=DARSIN (D) [D|>1 DA=0 D 
DA=DARCOS (D) 
268 DA=DTAN (D) |[D|22°*r DA=1 D 
DA=DCOTAN (D) 
269 DA=DTAN (D) D is too close to an DA=* D 
odd multiple of z 
DA=DCOTAN (D) D is too close to a DA=* D 
multiple of x 


For errors 271-275, let C=Xi+iXe 
Z=CEXP (C) 


















Xi > 174.673 Z=* (COS X:+SIN Xz) 






Z=CEXP (C) |X2|=2'*ar Z=0+0i 
Z=CLOG (C) C=0+4+0i Z=—*+0i 
Z=CSIN (C) \X|=2* 9 Z=0+0i 
Z=CCOS (C) 










Z=CSIN (C) Xo> 174.673 Z=*(SIN X:+iCOS X:) 






















2 
Z=CCOS (C) Z=*(COS X.-iSIN X,) 
2 
Z=CSIN (C) X2<—174.673 Z=*(SIN X.—iCOS X:) 
2, 










Z=CCOS (C) =*(COS Xi+ iSIN Xi) 


2 




















For errors 281-285, let CD =Xi+ixX:2 
Z=CDEXP (CD) 










Xi >174.673 Z=*(COS X2+iSIN Xe) 









Z=CDEXP (CD) |X2|=2°%r Z=0+0i 
Z=CDLOG (CD) CD=0+0i Z=—*+0i 
Z=CDSIN (CD) [Xs] 2° Z=0+0i 


Z=CDCOS (CD) 
Z=CDSIN (CD) 









X:>174.673 Z=*(SIN X,+iCOS X:) 
2 

Z=* (COS X:—iSIN X,) 
2 

Z=*(SIN X:—iCOS X:) 
2 

Z=*(COS X, +iSIN X:) 










Z=CDCOS (CD) 











Z=CDSIN (CD) Xe< — 174.673 















Z=CDCOS (CD) 











Y=GAMMA (X) X22 or 















X2257.5744 
Y=ALGAMA (X) X<0 or 

X24.2937*10" 
DA=DGAMMA (D) Ds2°” or 

D2257.5774 
DA=DLGAMA (D) D=0 or 






D24.2937* 10” 






Type 
I, J Variables of INTEGER*4 

X, XA, Y Variables of REAL*4 

D, DA, DB Variables of REAL*8 

Cc, CA Variables of COMPLEX*8 

Z, Xi, Xe Complex variables to be given the length of the functioned argument when they appear 
CD Variables of COMPLEX*16 













1. The user-supplied answer is obtained by recomputation of the function using the value set by the user routine for the parameters 
listed. (See “User-Supplied Corrective Action.” ) 
2. The largest number that can be represented in floating point is indicated by *. 
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Table 12. Corrective Action After Program Interrupt Occurrence 


Program Interrupt Messages 


Reason for 
Interrupt 
(See Note 1) 


Exponent overflow 
(Interrupt Code 12) 


Parameters 
Passed to 
User Exit 


D,I 


None 


Error 
a 


tee underflow 
(Interrupt Code 13) 

Divide check, 

Integer divide 

(Interrupt Code 9), 

Decimal divide 

(Interrupt Code 11), 

Floating point divide 
(Interrupt Code 15). 

See Note 4. 

Specification interrupt 

‘(Interrupt Code 6) occurs 

for boundary misalignment. 

Other interrupts occur 

during boundary alignment 

adjustment. They will be 

shown with this error 

code, and the PSW portion 

of the message will 

identify the interrupt. 


Type 
A variable REAL*8 





Standard User-Supplied 
Corrective Corrective 
Action Action 


User may alter D. 
See Note 2. 


Result register set to 
the largest possible 
floating-point number. 
The sign of the result 
register is not altered. 


The result register is User may alter D. 
set to zero. See Note 2. 


There is no standard See Note 3. 
fixup. Result registers 

No special corrective See Note 3. 
action other than correcting 

boundary misalignments. 


are not touched. 


Description 


This variable contains the contents of the result register after 
the interrupt. 

This variable contains the “exponent” as an integer value for the 
number in D. It may be used to determine the amount of 
the underflow or overflow. The value in I is not the true expo- 
nent, but what was left in the exponent field of a floating-point 


A variable INTEGER*4 





number after the interrupt. 


Notes to Table 12 


Note 1: A program interrupt occurs asynchronously. 
Interrupts are described in IBM System/360 Operating 
System: Principles of Operation, Form A22-6821. 


Note 2: The user-exit routine may supply an alternate 
answer for the setting of the result register. This is 
accomplished by placing a value for D in the user-exit 
routine. Although the interrupt may be caused by a 
long or short floating-point operation, the user-exit 
routine need not be concerned with this. The user-exit 
routine should always set a REAL*8 variable, and the 
FORTRAN library will load short or long depending 
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upon the floating-point operation that caused the in- 
terrupt. 


Note 3: The user-exit routine does not have the abil- 
ity to change result registers after a divide check. The 
boundary alignment adjustments are informational 
messages and there is nothing to alter before execution 
continues. 


Note 4: For floating-point divide check, the contents 
of the result register is shown in the message. 


The service subprograms supplied in the FORTRAN 
library are divided into two groups: one group tests 
machine indicators and the other group performs 
utility functions. Service subprograms are called by 
using the appropriate entry name in a FORTRAN source 
language CALL statement. 


Machine Indicator Test Subprograms 


The machine indicator subprograms (xxxFsLIT, 
XXXFOVER, and xxxFpvcH) test the status of pseudo in- 
dicators and may return a value to the calling program. 
When the indicator is zero, it is off; when the indica- 
tor is other than zero, it is on. In the following descrip- 
tions of the subprograms, i represents an integer 
expression and j represents an integer variable. 


xxxFSLIT Subprogram 


The xxx¥sLir subprogram is used to alter, test, and/or 
record the status of pseudo sense lights. Either of two 
entry names (SLITE or SLITET) is used to call the sub- 
program. The particular entry name used in the CALL 
statement depends upon the operation to be per- 
formed. 

If the four sense lights are to be turned oFF or one 
sense light is to be turned on, entry name SLITE is used. 
The source language statement is: 


CALL SLITE(?) 


where i has a value of 0, 1, 2, 3, or 4. 

If the value of i is 0, the four sense lights are turned 
off; if the value of i is 1, 2, 3, or 4, the corresponding 
sense light is turned on. If the value of i is not 0, 1, 2, 
3, or 4, an error message is issued and execution of 
this module (or phase) is terminated. (This error 
message is explained in Appendix C. ) 

If a sense light is to be tested and its status recorded, 
entry name SLITET is used. After a sense light is tested, 
it is set to off. The source language statement is: 


CALL SLITET (i, 7) 
where: 


i has a value of 1, 2, 3, or 4, and indicates which 
sense light to test. 

j is set to 1 if the sense light was on; or to 2 if the 
sense light was off. 


If the value of i is not 1, 2, 3, or 4, an error message 
is issued and execution of this module (or phase) is 


terminated. (This error message is explained in Appen- 
dix C.) 


Service Subprograms 


xxxFOVER Subprogram 
The xxxFover subprogram tests for an exponent over- 
flow or underflow exception and returns a value that 
indicates the existing condition. After testing, the over- 
flow indicator is turned off. This subprogram is called 
by using the entry name OVERFL in a CALL statement. 
The source language statement is: 
CALL OVERFL (/) 
where: 
j is set to 1 if a floating-point overflow condition 
exists; to 2 if no overflow or underflow condition 
exists; or to 3 if a floating-point underflow condi- 
tion exists. A detailed description of each exception 
is given in Appendix C. 


xxxFDVCH Subprogram 

The xxxFpvcH subprogram tests for a divide-check 
exception and returns a value that indicates the exist- 
ing condition. After testing, the divide-check indicator 
is turned off. This subprogram is called by using entry 
name DVCHK in a CALL statement. The source language 
statement is: 


CALL DVCHK (/) 
where: 


j is set to 1 if the divide-check indicator was on; or 
to 2 if the indicator was off. A detailed description of 
the divide-check exception is given in Appendix C. 


Utility Subprograms 

The utility subprograms perform two operations for 
the FORTRAN programmer: they either terminate execu- 
tion (xxxFExir) or dump a specified area of storage 
(xxxFDUMP ). 


xxxFEXIT Subprogram 


The xxxFexir subprogram terminates execution of this 
load module (or phase) and returns control to the 
operating system. (This subprogram performs a func- 
tion similar to that performed by the stop statement. ) 
The xxx¥exir subprogram is called by using the entry 
name ExIT in a CALL statement. The source language 
statement is: 


CALL EXIT 


xxxFDUMP Subprogram 


The xxxFpuMp subprogram dumps a specified area of 
storage. Either of two entry names (DUMP or PDUMP) 
can be used to call the subprogram. The entry name 
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is followed by the limits of the area to be dumped and 
the format specification. The entry name used in the 
CALL statement depends upon the nature of the dump 
to be taken. 

If execution of this load module (or phase) is to be 
terminated after the dump is taken, entry name DuMP 
is used. The source language statement is: 


CALL DUMP (a, bi, fi e+e Qn, ba, fn) 
where: 

a and b are variables that indicate the limits of stor- 
age to be dumped (either a or b may represent the 
upper or lower limits of storage). 

f indicates the dump format and may be one of the 
integers given in Table 13. The formats available 
depend upon the compiler in use. A sample printout 
for each format is given in Appendix F. 


Table 13. The xxxFDUMP Subprogram Format Specifications 


FORTRAN IV (E) FORTRAN IV 


0 specifies hexadecimal 0 specifies hexademical 
4 specifies INTEGER 1 specifies LOGICAL *1 
5 specifies REAL 2 specifies LOGICAL *4 
6 specifies DOUBLE 3 specifies INTEGER *2 
PRECISION 4 specifies INTEGER *4 

5 

6 

7 

8 




















specifies REAL *4 
specifies REAL *8 
specifies COMPLEX *8 
specifies COMPLEX *16 


9 specifies literal 


If execution is to be resumed after the dump is taken, 
entry name ppuMpP is used. The source language state- 
ment is: 


CALL PDUMP (a, bi, fi,..., dn, Dn, fn) 


where a, b, and f have the same meaning as explained 
previously. 


Programming Considerations 


A load module (or, in the Model 44 Programming 
System, a member of the phase library) may occupy 
a different area of storage each time it is executed. To 
ensure that the appropriate areas of storage are dumped, 
the following conventions should be observed. 


Note: In the following examples, a is a variable in 
COMMON, B is a real number, and the array TABLE is 
dimensioned as: 

DIMENSION TABLE (20) 


If an array and a variable are to be. dumped at the 
same time, a separate set of arguments should be 
used for the array and for the variable. The specifica- 
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tion of limits for the array should be from the first 
element in the array to the last element. For example, 
the following call to the mcrpump subprogram could 
be used to dump TABLE and B in hexadecimal format 
and terminate execution after the dump is taken: 


CALL DUMP (TABLE (1), TABLE (20), 0, B, B, 0) 


If an area of storage in comMON is to be dumped 
at the same time as an area of storage not in COMMON, 
the arguments for the area in CoMMON should be given 
separately. For example, the following call to the 
XXXFDUMP subprogram could be used to dump the 
variables a and B in REAL *8 format without terminat- 
ing execution: 

CALL PDUMP (A,A,6,B,B,6) 


If variables not in comMON are to be dumped, each 
variable must be listed separately in the argument 
list. For example, if R, p, and Q are defined implicitly 
in the program, the statement 


CALL PDUMP (R,R,5,P,P,5,0,0,5) 


should be used to dump the three variables. If the 
statement 
CALL PDUMP (R,Q,3) 


is used, all main storage between R and Q is dumped, 
which may or may not include p, and may include 
other variables. 

If an array and a variable are passed to a subroutine 
as arguments, the arguments in the call to the xxxFDUMP 
subprogram in the subroutine should specify the 


parameters used in the definition of the subroutine. 


For example, if the subroutine susi is defined as: 


SUBROUTINE SUBI (X, Y) 
DIMENSION X(10) 


and the call to susr within the source module is: 
DIMENSION A(10) 


CALL SUBI (A, B) 


then the following statement in the subroutine should 
be used to dump the variables in hexadecimal format 
without terminating execution: 


CALL PDUMP (X(1), X(10), 0, Y, Y, 0) 
If the statement 
CALL PDUMP (X(1), Y, 0) 


is used, all storage between a(1) and y is dumped, 
due to the method of transmitting arguments. 


Appendix A. Algorithms 


Appendix A contains information about the computations used in the explicitly 
called mathematical subprograms. This information is arranged in alphabetical 
order, according to the module (or phase) name of the subprogram. The entry 
names associated with each subprogram are given in parentheses after the module 
(or phase) name. 


The information for each subprogram is divided into two parts. The first part 
describes the algorithm used; the second part describes the effect of an argument 
error upon the accuracy of the answer returned. 


The presentation of each algorithm is divided into its major computational steps; 
the formulas necessary for each step are supplied. Some of the formulas are 
widely known; those that are not so widely known are derived from more common 
formulas. The process leading from the common formula to the computational 
formula is sketched in enough detail so that the derivation can be reconstructed 
by any one who has an understanding of higher mathematics and access to the 
common texts on numerical analysis. 


The accuracy of an answer produced by these algorithms is influenced by two 
factors: the performance of the subprogram (see Appendix B) and the accuracy 
of the argument. The effect of an argument error upon the accuracy of an answer 
depends solely upon the mathematical function involved and not upon the partic- 
ular coding used in the subprogram. 


A guide to the propagational effect of argument errors is provided because 
argument errors always influence the accuracy of answers whether the errors are 
accumulated prior to use of the subprogram or introduced by newly converted 
data. This guide (expressed as a simple formula where possible) is intended to 
assist users in assessing the effect of an argument error. 


The following symbols are used in this appendix to describe the effect of an 
argument error upon the accuracy of the answer: 








SYMBOL EXPLANATION 
q(x) The result given by the subprogram. 
f(x) The correct result. 
: f(x) — g(x) The relative error of the result 
f(x) given by the subprogram. 
8 The relative error of the argument. 
E | f(x) — g(x) | The absolute error of the result 
given by the subprogram. 
A The absolute error of the argument. 


The notation used for the continued fractions complies with the specifications 
set by the National Bureau of Standards.? 


1Any of the common numerical analysis texts may be used as a reference. One such text is F. Hildebrand’s 
Introduction to Numerical Analysis (McGraw-Hill Book Company, Inc., New York, N, Y., 1956). Background 
information for algorithms that use continued fractions may be found in H. S. Wall’s Analytic Theory of 
Continued Fractions (D. VanNostrand Co., Inc., Princeton, N. J., 1948). 

2¥For more information, see Milton ‘Abramowits and Irene A. Stegun (editors), Handbook of Moriemation! 
Functions, Applied Mathematics Series-55 (National Bureau of Standards, Washington, D.C.), 1965. 


Appendix A. Algorithms 21 


xxxCLABS (CDABS) and xxxCSABS (CABS) Subprograms 
1. Write |x + iy] = a + ib. 
2. Ifx = y = 0, thena = 0 and b = 0. 


3. Let vo, = max (|x|, |y|), and 
v2 = min (|x|, |y]). 


[fay VS 
Then aie 1+ (=) ,and b = 0. 
: 1 








? 


The algorithms for both complex absolute value subprograms are identical. 
Each subprogram uses the appropriate real square root subprogram (xxxLSQRT 
Or XXXSSQRT ). 


xxxCLEXP (CDEXP) and xxxCSEXP (CEXP) Subprograms 


Algorithm 


The value of e”+” is computed as e” * cos(y) + i* e” * sin(y). The algorithms 
for both complex exponential subprograms are identical. Each subprogram uses 
the appropriate real exponential subprogram (xxxLExP or xxxsExP ) and the appro- 
priate real sine/ cosine subprogram (xxxLSCN or XXXSSCN). 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If e*+” = R » e, then H = y and e(R) ~ A (x). 


xxxCLLOG (CDLOG) and xxxCSLOG (CLOG) Subprograms 


Algorithm 
1. Write log. (x + iy) = a + ib. 


2. Then, a = log, |x + iy|and b = the principle value of arctan 2 : 


The algorithms for both complex natural logarithm subprograms are identical. 
Each subprogram uses the appropriate complex absolute value subprogram 
(XxxCLABS or xxxcsABs), the appropriate real natural logarithm subprogram 
(XxxLLOG or xxxsLoGc), and the appropriate arctangent subprogram (xxxLATN2 
OY XXXSATN2 ). 


Effect of an Argument Error 


The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = r* e” and log, (x + iy) = a + ib, then h = b 
and E(a) = 8(r). 


xxxCLSQT (CDSQRT) and xxxCSSQT (CSQRT) Subprograms 
Algorithm 

l. Write Ve + iy =a + ib. 

2. Ifx = y = 0, thena = 0Oandb = 0. 


3. If 2 0,thena = | + |x + iy 
- 2 
and b = +, 

2a 


22, 


4. Ifx < 0, then b = (sign y). (te 
da=/. 
and ¢ = = 


The algorithms for both complex square root subprograms are identical. Each 
subprogram uses the appropriate real square root subprogram (xxxLsQrT or 
XXXSSQRT ). 


Effect of an Argument Error 
The effect of an argument error depends upon the accuracy of the individual 
parts of the argument. If x + iy = re e” and Vx + iy = R° e'*, then 


Re = a(r), and e( H) ~ 8(h). 


xxxCLSCN Subprogram (CDSIN and CDCOS) 
Algorithm 
]. If the sine is desired, then 

sin(x + iy) = sin(x) * cosh(y) + i* cos(x) ° sinh(y). 

If the cosine is desired, then 

cos(x + iy) = cos(x) * cosh(y) — i* sin(x) * sinh(y). 

2. Ifx < 0, then sinh( —x) = —sinh(x). 
er —-— 


3. If x > 0.3465736, then sinh(x) = a 





4. If0 = x S 0.3465736, then compute sinh (x) by use of the polynomial: 
sins) = Ay + yx? + agxt +... + asx. 

The coefficients are obtained by expanding the polynomial with respect to the 

Chebyshev polynomials over the range 0 = x? = 0.120113. The relative error 

of this approximation is less than 2~21-8, 


5. The value of cosh(x) is computed as cosh(x) = sinh|x| + i ‘ 
This computation uses the real exponential subprogram (xxxLExp) and the 


real sine/ cosine subprogram (xxxLSCN ). 


Effect of an Argument Error 


To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument error in the xxxLscn, 
XXXLEXP, and xxxLSNH subprograms. 


xxxCSSCN Subprogram (CSIN and CCOS) 


Algorithm 
1. If the sine is desired, then 
sin (x + iy) = sin(x) * cosh(y) + i* cos(x) * sinh(y). 
If the cosine is desired, then 
cos(x + iy) = cos(x) * cosh(y) — i* sin(x) * sinh(y). 
2. Ifx < 0, then sinh( —x) = —sinh(x). 


Ee? —-— 


3. Ifx > 0.3465736, then sinh(x) = 5 a ; 
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4. If0 =x = 0.3465736, then compute sinh (x) by use of the polynomial: 
sae = Ay + ax? + dox'. 

The coefficients are obtained by expanding the polynomial with respect to the 

Chebyshev polynomials over the range 0 S x? < 0.120113. The relative error 

of this approximation is less than 2-7-4, 


5. The value of cosh(x) is computed as cosh(x) = sinh|x| + = ; 
ele 
This computation uses the real exponential subprogram (xxxsExp) and the 
real sine/cosine subprogram (xxxsscn ). 


Effect of an Argument Error 


To understand the effect of an argument error upon the accuracy of the answer, 
the programmer must understand the effect of an argument in the xxxsscn, XXxxSEXP, 
and xxxsscnH subprograms. 


xxxLASCN Subprogram (DARSIN and DARCOS) 


Algorithm 
1lfOscs + , then compute arccos(x) as: 
arccos(x) = or — arcsin(x). 
Ifosxcs + , then compute arcsin(x) by a polynomial of the form: 
arcsin(x) = x + ex? + cox® +... A C12x®. 
; : ; ; arcsin(x) 
The coefficients are obtained by expanding the function f(z) = eee 
1 


z = x*, with respect to the Chebyshev polynomials over the range, 0 = x ST: 


The relative error of this approximation is less than 2—55-7, 
2. If +. < x <1, then compute arcsin (x) as: 


arcsin(x) = —-— arccos(x). 


pi 
2 
If +. < x <1, then compute arccos (x) as: 


[j—-* 
arccos(x) = 2° arcsin Cy : 3 *). 


This case is now reduced to the first case because within these limits, 





ee 


9S ya Sr 


This computation uses the real square root subprogram (xxxLSQRT ). 
3. If -1S-x < 0, then aresin(x) = —arcsin|s| 
and arccos(x) = 7 — arccos|x]. 


This reduces these cases to one of the two positive cases. 


Effect of an Argument Error 
tA 
Vike x 
a small A causes a substantial error in the answer. For the arcsine, « ~ 8 if the 
value of x is small. 


EK ~ 





. For small values of x, E ~ A. Toward the limits (+1) of the range 


IHCLATAN Subprogram (DATAN) 


Algorithm 
1. Reduce the computation of arctan (x) to the caseO0 Sx <1 by using 
arctan( —x) = —arctan(x) or 
1 7 
——— TT —_—_ — t 
arctan ry 5 arctan |x| 


2. If necessary, reduce the computation further to the case |x| = tan 15° by using 
arctan (x) = 30° + arctan ee 
x+ V3 


The value of Ri Sea <= tan 15° if the value of x is within the range, 








tan 15° <x<1. The value of (VW3*x —1) is computed as 


] 


(V3 —-1)x- 7 -+ + x to avoid the loss of significant digits. 


3. For |x| S tan 15°, use a continued fraction of the form: 


arctan (x) es a,x? As as a4 
x _ (by + x?) + (bg +x?) + (bg + x?) + (bg t+ x?) +07" 


The relative error of this approximation is less than 2-57. The coefficients of 
this formula were derived by transforming the continued fraction: 


i 304 16 * 25 
arctan (x) _ 1+ 3 25° 7 7° 8l¢11 
x . ($+=)- (s<9% 3") i (saat 4) 
5 5°9 9°13 
4*3°9 
5°11 * 169. 


3°37 a 
(Gssa7 +? ) - 
where w has an approximate value of x —? + 40) but the true 


64° 27 
5 * 289 * 19 


value of w is vee 
(Gente) 


9 
5°ll- 13-17 





3°7°17 


Effect of an Argument Error 


a a 
Lx 
effect of « upon § diminishes. 


For small values of x, « ~ 8, and as the value of x increases, the 


xxxLATN2 Subprogram (DATAN and DATAN2) 
Algorithm 


Tv 


7° 





aoe 


2 


1. For arctan(%;, x2), if either x, = 0 or > 25° the answer = (sign x) ° 








Otherwise, if x. > 0, the answer = arctan (= ) , and 
2 


if x2 < 0, the answer = arctan (2) + (sign x,) ° 7. 
2 
The rest of the computation is identical for either one or two arguments. 
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2. Reduce the computation of arctan(x) to the case OS x S 1, by using 


arctan( —x) = —arctan(x), or 
arctan oe Ve arctan|+x| 
Ix] 72 


3. If necessary, reduce the computation further to the case |x| = tan 15° by using 


arctan(x) = 30° + arctan vee 
x+ V3 
V3°x-—1 
x+ V3 
tan 15° < x <1. The value of (\/3* x — 1) is computed as (3 — 1) x - 1 
to avoid the loss of significant digits. 


The value of < tan 15° if the value x is within the range, 








4. For |x| S tan 15°, use a continued fraction of the form: 
arctan(x) oe ax? do a3 a4 ae 
* (by + x?) + (bp +x?) + (bg +x?) + (ba +x?) + 
The relative error of this approximation is less than 2—57-°, The coefficients of 
this formula were derived by transforming the continued fraction: 


=1 34 16 + 25 
arctan(x) esse 3 25° 7 7°8le ll 
x 3 ee Nes 23 =. ae 59 > 
(s+ ) (sat? ) (seat? 
4°3°9 
5° 11° 169 





i ay 
(get? ) i 


2 —— 
where w has an approximate value of Sell? 13°17 (—x—2 + 40) but the 
64 * 27 
5° 289+ 19 


] is ————_——_—__ .. . 
true value of w is i79 ae 
“4 3¢7¢17 
Effect of an Argument Error 


E~ . For small values of x, « ~ 8, and as the value of x increases, the effect 


see. 
LP x? 
of e upon § diminishes. 


xxxLERF Subprogram (DERF and DERFC) 





Algorithm 
1. If0 Sx < I, then compute the error function by the following approximation: 
erf(x) s2 x (do + ayx? + doxt +... + yx”). 
The coefficients were obtained by expanding the function f(z) = wil o). 


x 
z = x?, with respect to the Chebyshev polynomials over the range, 0 = x < 1. 


The relative error of this approximation is less than 1.07 * 2-5". The value of 
the complemented error function is computed as erfc(x) = 1 — erf(x) and is 


1 
greater than 16° 


2. If l= x < 2.0400009, then compute the complemented error function by the 
following approximation: 


erfc (x) > bo + biz + bez? +...+ byigz18 
where z = x — To and Ty = 1.999999. The coefficients were obtained by ex- 


panding the function f(z) = erfc(z + To) with respect to the Chebyshev poly- 
nomials over the range —1 = x = 0.04. The absolute error of this approxima- 
tion is less than 1.5 * 2-6. The limits of this range and the base value for To 

~were used to minimize the hexadecimal truncation error. The value of the com- 
plemented error function within this range is greater than = . The value of 
the error function is computed as erf(x) = 1 — erfe(<x). 


3. If 2.0400009 = x = 13.306, then compute the complemented error function by 

the following approximation: 
6S. (Co + C1x-? + cox* +... + Caox—*) C-® 

x 

The coefficients were obtained by expanding the function 
f(z) = erfc (x) *x* e”, z = x—%, with respect to the Chebyshev polynomials 
over the range 2.04-2 > z = 13.306—?. The relative error of this approximation 
ranges from 2-53 at 2.04 to 2-5! at 13.306. This computation uses the real ex- 
ponential subprogram (xxxLExP ). 
If x <= 6.092, then the error function is computed as erf(x) = 1 — erfc(x). 
If x > 6.092, then the error function is = 1. 


4, If 13.306 < x, then the error function is = 1, and the complemented error 
function is = 0. 


erfc (x) 


5. If x < 0, then reduce to a case involving a positive argument by the use of 
the following formulas: 
erf(—x) = —erf(x) and erfe(—x) = 2 — erfe(x). 


Effect of an Argument Error 


E ~ e-” A. For the error function, as the magnitude of the argument exceeds 
1, the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, « ~ 8. For the complemented error function, if the value of 


—g? 


x is greater than 1, erfc(x) ~ x=: Therefore, « ~ 2x° ¢ 8. If the value of x 
x 


is negative or less than 1, then e ~ e~”* A, 


xxxLEXP Subprogram (DEXP) 

Algorithm 

1. If x < — 180. 2183, then 0 is given as the answer. 
2. Divide x by log.2 and write 


x c 
= = (4a —- b -——-—d 
Y= Tea | ie” 

where a, b, and c are integers, O= b<3, OSc=15, and d is within the range 


O<de = Then, e? = 2 = 162+9->« 9-c/16. 9-4 





3. Compute 2—¢ by using the Chebyshev interpolation of degree 6 over the range, 
0Osd< = . The maximum relative error of this computation is 2-57. 


4. Ifc > 0, then multiply 2-4 by 2-°/18, (The 15 values of 2-°/16 for lS c 15 
are included in the subprogram. ) 


5. If b > 0, then halve the result b times. 
6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 


Effect of an Argument Error 


E ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8° x. 
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~ xxxLGAMA Subprogram (DGAMMA and DLGAMA) 


Algorithm 

1. 10 < x S 2-252, then compute log-gamma as log.'(x) = —log.(x). This 
computation uses the real logarithm subprogram (xxxLLoc). 

2. If 2-252 < x < 8, then compute log-gamma by taking the natural logarithm 
of the value obtained for gamma. The computation of gamma depends upon 
the range into which the argument falls. 


3. If 2-252 < x < I, then use (x) = a to reduce to the next case. 


4. If 1 = x S 2, then compute gamma by the following approximation: 
T(x) = do + 12 + ez? +... + da92” 

where z = x — 1.5. The coefficients were obtained by expanding the function 
f(z) = I(x) with respect to the Chebyshev polynomial for |z| S 0.5. The 
absolute error of this approximation is less than 1.5 + 2758, 

5. If2 <x < 8, then use I(x) = (x — 1) I(x — 1) to reduce to the preced- 
ing case. 

6. If 8 = x, then compute log-gamma by the use of Stirling’s formula: 


ibaa) Se aloeay = 1).= $ loge(2) 5 log.(2r) + G(x). 


The modifier term G(x) is computed as 
G(x) = Dix? a box 8 - byx—5 “+ byx-7 oi bsx-°®. 





The coefficients were obtained by expanding the function f(z) = ots) ; 
Z = x-*, with respect to the Chebyshev polynomials over the range 


0 < z < 8-?. The absolute error of the approximation for G(x) is less than 
x * 2-58 Because, in this range, x < log.I'(x), the contribution of this error 
to the relative error of the value for log-gamma is less than 2—5*. This com- 
putation uses the real logarithm subprogram (xxxLLoc). 

For gamma, compute I'(x) = e’, where y is the value obtained for log- 
gamma. This computation uses the real exponential subprogram (xxxLEXP). 


Effect of an Argument Error 
e ~ ¥(x) * A for gamma, and E ~ y(x) * A for log-gamma, where y is the 
digamma function. 

If 3 <x < 3, then —2 < ¥(x) < 1. Therefore, E ~ A for log-gamma. 


However, because x = 1 and x = 2 are zeros of the log-gamma function, even 
a small § can cause a substantial « in this range. 

If the value of x is large, then y(x) ~ log.(x). Therefore, for gamma, 
e~ 8*x° log.(x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x,« ~ 8. 


xxxLLOG Subprogram (DLOG and DLOG10) 


Algorithm 
1, Write x = 16? + 2—%* m, where p is the exponent, q is an integer, 0S q $3, 
and m is within the range, = =m<l. 
2. Define two constants, a and b (where.a = base point and 2~ = a) as follows: 
If + =sm< then a = Sand b =], 
If en <1,thena = landb = 0. 


+z 


Leen - and |z| < 0.1716. 


mta 
4, Now, x = 24?-9-2 -+ 


5. Finally, log, (; = =) is computed by using the Chebyshev interpolation of 


degree 7 in z? over the range, 0 z? 0. 02944. The maximum relative error 
of this approximation is 2—5°-6, 








. Then, m = a* = 


and log.x = (4p — q — b) log.2 + los. ( 


3. Write z = 


+ 2 
— 2’ 


1+2z 
1—z/s° 











6. If the common logarithm is desired, then logigx = logice ° log.x. 


Effect of an Argument Error 


E ~ 8. Therefore, if the value of the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 


xxxLSCN Subprogram (DSIN and DCOS) 
Algorithm 


1. Divide |x| by a and separate the quotient (z) into its integer part (q) and 





4 ; 
its fraction part (r), Then, z = |x = q + r, where q is an integer and r 
TT 





is within the range, OS7 <1. 


2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of gq reduces the general case to the computation of 
sin (x) for x 2 0, because 


cos (+x) = sin (i: +5) , and 


sin (—x) = sin (|x| + 7). 
3. Let go = q mod 8, 


Then, for qo = 0, sin (x) = sin G r) 


Fe eee (4a = *)) 
rr) 

qo = 3, sin (x) = sin (qa — ")) 
qo = 4,sin (x) = —sin (+"*) 
pes win (4 ene (qa = r)) 
i eG ane Sele (+ r) 
goin eae (+0 = r)) 


These formulas reduce each case to the computation of either sin (+ r) or 





qo = 2, sin (x) = cos ( 


4 


7 * 11 }; where 1, is either r or (1 —r), and is within the range, OSS. 





COs 





zi * 1, } or cos aT r; } is computed, using the Chebyshev 


interpolation of degree 6 in rj? for the sine, and of degree 7 in 1,? for the cosine. 
The maximum relative error of the sine polynomial is 2—5* and that of the 
cosine polynomial is 2—&4:8, _ 


4, Finally, either sin 
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Effect of an Argument Error 


E ~ A. As the value of the argument increases, A increases. Because the function 
value diminishes periodically, no consistent relative error control can be main- 


tained outside of the principal range, — o. oe Re os ' 


xxxLSCNH Subprogram (DSINH and DCOSH) 


Algorithm 
1. If |x| < 0.3465736, then compute sinh(x) as: 
sinh(x) = x + cx? + cox® + cgx? + cgx® + coxh. 


The coefficients are obtained by expanding the function f(z) = = , 


% = x?, with respect to the Chebyshev polynomials over the range, 
0=z < 0.12011326. The relative error of this approximation is less than 2—§-9, 
2. If either |x| = 0.3465736 or the cosh(x) is desired, obtain w = e'!, Then, 
—1 panes —1 
oe and sinh(x) = (sign x) eee The real expo- 


nential subprogram (xxxLExP ) is used to compute the value of w. 


cosh(x) = 


Effect of an Argument Error 
For the hyperbolic sine, E ~ A * cosh(x) and « ~ A * coth(x). 
For the hyperbolic cosine, E ~ A * sinh(x) and « ~ A * tanh(x). 
Specifically, for the cosine, E ~ A over the entire range; for the sine, «e ~ 8 
for the small values of x. 


xxxLSQRT Subprogram (DSQRT) 


Algorithm 


1. If x = 0, then the answer is 0. 


2. Write x = 16??—%* m, where 2p — q is the exponent and q equals either 0 or 1; 
m is the mantissa and is within the range, 4 <m<l. 

3, Then, Vx = 167+ 2-24+ W/m. 

4. For the first approximation of \/x , compute the following: 


Yo = 2-74 16? + (s+ 2 ‘m). 





9 


The maximum relative error of this approximation is 9° 


5. Apply the Newton-Raphson iteration 


1 bs 
Bp rahe) 


four times to Yo (the first two times in the short form and the last two times in 
the long form). The final step is performed as 


1fx 
— + —j— - 
Ys = Ys 5) ( Ys us) 
to minimize the computational truncation error. The maximum relative error 
of the final result is theoretically 2-79, 


Effect of an Argument Error 


e~ —§ 


xxxLTANH Subprogram (DTANH) 


Algorithm 


1. 


2. 


3. 
4, 


If |x| < 0.54931, then use the following fractional approximation: 


tanh (x) L- Ax? + dox* + agx® + x8 
Xx bo + byx? + boxt + bax® + x8 
where: 
a, = 676440.765 bo = 2029322.295 
a3 = 594.459 b. = 52028.55 


b; = 630.476 


The maximum relative error of this approximation is 2~*5, The formula was 
obtained by transforming the continued fraction 





tanh (x) _ 1+ , ea i 
x 3+ 5+°°° 15+w 
where w has an approximate value of 0.017, but the true value of w is 
, ane 
17+ 19+ °° 


If 0.54931 = x < 20.101, then use the identity tanh (x) = 1 — This 


e* + 1° 
computation uses the double precision exponential subprogram (xxxLEXP). 
If x = 20.101, then tanh (x) = 1. 


If x = — 0.54931, then use the identity tanh (x) = — tanh ( —x). 


Effect of an Argument Error 


E 


~ (1 — tanh? x) A, ande ~ 


2A 


sigh (20) For small values of x , « ~ 8. As the 


value of x increases, the effect of § upon « diminishes. 


xxxLTNCT Subprogram (DTAN and DCOTAN) 


Algorithm 


I. 


. The values of tan Ze w) and cot ( 7 


Divide |x| by and separate the result into the integer part (q) and the 


fraction part (r). Then, |x| = z (q +f). 


. Obtain the reduced argument (w) as follows: 


if q is even, then w = r. 
if q is odd, then w = 1 — r. 
The range of the reduced argument is 0 S w S 1. 


. Let do = q mod 4. 








Then, for go = 0, tan |x| = tan ( a “Ww ) and cot |x| = cot (4 ° w) 
<0 
4 


qo = 1, tan |x| = cot G -w ) and cot |x| = tan ( ‘w) 
do = 2, tan |x| = —cot ( = -w) and cot |x| — tan ( 








Fags 
4 4 


go = 3, tan [x| = tan ( 7 -w) and cot |x| = —cot ( i -w) 














( a ° w) are computed as the ratio 


of two polynomials. 
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7 _ we P(w?) : — __O(w*) _ 
tan (4 w) = On)? and cot (5 w) = w+ P(w?) 
where P(w?) is of degree 3 and Q(w?) is of degree 4 in w®. The etic 
of P and Q are obtained by economizing the continued fraction 





tan(z) _ 1- ze gt ge 
z Gan Foe, Geo 
in the following way. 
a ae ae 2 2 ‘ 
Write. @2@) 2,9 # # 2@ # z z z 


38—- 5-— 7-— 9- (11+ d,)— (13+ d,.)— (15 4+ ds) 
and determine the values for d,, dz, and ds so that the right-hand expression 
gives the exact answers for z? = 0.395, 0.542, and 0.607. Then the maximum 


relative error of this formula over the range 0 S z S— is 3.4 + 10-19, 


4 
Change the variable from z to w =7 * z and rewrite the formula to obtain 
P(w*) and QO(w?). 
3. If x < 0, then tan(x) = —tan |x|, and cot(x) = —cot [x]. 


Effect of an Argument Error 


, and e ~ for tan(x). Therefore, near the singularities 


Se 
cos? (x) 


2 
sin (2x) 
] ae: 
of x = ( k + 9 I where k is an integer, no error control can be maintained. 


This is also true for cotan(x) for values of x near kz, where k is an integer. 


xxxSASCN Subprogram (ARSIN and ARCOS) 
Algorithm 


Lifos*es >, then compute arccos(x) as: 


arccos(x) = el arcsin(x). 


If0sxcs > then compute arcsin(x) by a polynomial of the form: 


arcsin(x) 22x + yx? + cox® + cgx7 + c4x® + 5x11. 
: , : : arcsin (x) 
The coefficients are obtained by expanding the function f(z) = ae 


: ‘ 1 
% = x?, with respect to the Chebyshev polynomials over the range 0 = z =T° 
The relative error of this approximation is less than 2-27-5, 


2. If >< x <= I, then compute arcsin(x) as: 


arcsin(x) = o — arccos(x). 


If + < x S1, then compute arccos(x) as: 





arccos(x) = 2° arcsin : > *). 


This case is now reduced to the first case because within these limits, 





This computation uses the real square root subprogram (xxxssQrr ). 


3. If -1Sx < 0, then arcsin(x) = — arcsin |x|, and arccos(x) = a — arccos |x]. 


This reduces these cases to one of the two positive cases. 


Effect of an Argument Error 
A nme 

E ~ ee For small values of x, E ~ A. Toward the limits (+1) of the 
— 


range, a small A causes a substantial error in the answer. 


IHCSATAN Subprogram (ATAN) 
Algorithm 


1. Reduce the computation of arctan (x) to the case 0 = x S 1, by using 
arctan (—x) = —arctan (x), or 


arctan 2 ane arctan |x| 
|x| 2 


2. If necessary, reduce the computation further to the case |x| S tan 15° by using 


arctan (x) = 30° + arctan Sn) : 
at x+vV3 
V38ex-1 , es 
The value of |}~——~———} S tan 15° if the value of x is within the range, 
x+ V3 








tan 15° < x <1. The value of (V3 + x — 1) is computed as 
(V3 —1) x — 1+ x to avoid the loss of significant digits. 


3. For |x| S tan 15°, use the approximation formula: 


arctan (x) _ 0.55913709 
Se 0.60310579 — 0.05160454x? + 4+ 14087812 


This formula has a relative error less than 2~27-1 and can be obtained by 
transforming the continued fraction 


x2 
arctan (x) _ 1 = ee 5B 
x 38+ ( 5 -) 
rate ao — w 
t 
where w has an approximate value of ( 2 aii a 10-4, but the true 
4°5 
7*7°9 


value of w is 





f. Boye 
(fat? )+ 


The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 


Effect of an Argument Error 


Ew . For small values of x, « ~ 8; as the value of x increases, the effect 


2k 
Li x 
of § upon e diminishes. 


xxxSATN2 Subprogram (ATAN and ATAN2) 


Algorithm 
Xy 


2 


Otherwise, if x. > 0, the answer = arctan (=), and 
2 


> 924 the answer = (sign x1) * —. 


l. For arctan(x,, x2), if either x. = 0 or 5) 








if x. < 0, the answer = arctan (=) + (sign x;) 7. 
2 


The rest of the computation is identical for either one or two arguments, 
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2. Reduce the computation of arctan(x) to the case 0 = x = 1, by using 
arctan( —x) = — arctan(x) or 


TN x 
arctan (a)= 9 arctan |x|. 


3. If necessary, reduce the computation further to the case |x| S tan 15° by using 


arctan(x) = 30° + arctan Mose 
x+ V3 
V3ex-1 , ‘ oe 
The value of ea < tan 15° if the value of x is within the range, 
x 








tan 15° < x <1. The value of (V3 * x — 1) is computed as 
(V3 — 1) x — 1+ <x to avoid the loss of significant digits. 


4. For |x|= tan 15°, use the approximation formula: 


arctan(x) = js 0.55913709 
=a 0.60310579 — 0.0516045x +4 14087812 + 14087819. 


This formula has a relative error less than 2—27-! and can be obtained by trans- 
forming the continued fraction 


x? 
arctan(x) _ L- - et 
x 7 3+ ( 5 ~*) 
oor © — Ww 
7 
where w has an approximate value of ( _ Be + aT) 10-4 but the true 
4° 
value of w is EAS 
43 





7°11 
The original continued fraction can be obtained by transforming the Taylor 
series into continued fraction form. 


+ =*) + 


Effect of an Argument Error 


Ew~ . For small values of x, « ~ 8; as the value of x increases, the effect 


AT 
1+ x? 


-of § upon ¢ diminishes. 


xxxSERF Subprogram (ERF and ERFC) 
Algorithm 


1. If 0 S x = 1.317, then compute the error function by the following approxi- 
mation: 
erf(x) = x(do + a,x? + doxt +... + agx}?). 


The coefficients were obtained by expanding the function f(z) = ont) 2S, 


with respect to the Chebyshev polynomials over the range 0 S x = 1.317. The 
relative error of this approximation is less than 2-24. The value of the comple- 
mented error function is computed as erfc(x) = 1 — erf(x) and is greater 


l 
than 6° 


2. If 1.317 < x S 2.0400009, then compute the complemented error function by 
the following approximation: 
erfc(x) = bo + biz + bez? +... + D727 
where z = x — Ty and Ty) = 2.0400009. The coefficients were obtained by 
expanding the function f(z) = erfe(x + T ) with respect to the Chebyshev » 


polynomials over the range (1.317 — Ty) < z S 0. The absolute error of this 
approximation is less than 1.3 * 2—%°. The value of the complemented error 
function within the range 1.317 < x S Ty is greater than = . The value of the 
error function is computed as erf(x) = 1 — erfe(x). 


3. If To < x S 13.306, then compute the complemented error function by the fol- 

lowing approximation: 
(eg F ew=? + coxa +... F eee?) e- 

% 

The coefficients were obtained by expanding the function 
f(z) = erfc (x) * x e”, z = x—, with respect to the Chebyshev polynomials 
over the range T)~2 > z = 13.306—2. The relative error of this approximation 
is less than 1.2 + 2-23, This computation uses the real exponential subprogram 
( XxXSEXP ). 
If x = 3.9192, then the error function is computed as erf(x) = 1 — erfc(x). 
If x > 3.9192, then the error function is = 1. 


erfc(x) 


4. If 13.306 < x, then the error function is = 1, and the complemented error func- 
tion is = 0. 


5. If x < 0, then reduce to a case involving a positive argument by the use of the 
following formulas: 
erf(—x) = —erf(x) and 
erfc(—x) = 2 — erfe(x). 


Effect of an Argument Error 


E ~ e—” + A. For the error function, as the magnitude of the argument exceeds 1, 
the effect of an argument error upon the final accuracy diminishes rapidly. For 
small values of x, « ~ 8. For the complemented error function, if the value of x is 


greater than I, erfc(x) ~ — . Therefore, « ~ 2 x? « 8. If the value of x is negative 


or less than 1, then e ~ e~-” °° A. 


xxxSEXP Subprogram (EXP) 
Algorithm 
1. Ifx < —180.218, then 0 is given as the answer. 
2. If |x| < 2-8, then 1 is given as the answer. 
3. Otherwise, divide x by log.2 and write 
¥ = og = (4a — b — d) 
where a and b are integers,O = b=3and0Sd< 1. 
Then, e”? = 2Y = (167+ 2->+2-4), 
4, Compute 2~¢ by the following fractional approximation: 


gta — [ 2d 


817 97297 \ 
ee d*? + d + 9.9545958 — 


& + 87417497 J 


This formula can be obtained by transforming the Gaussian-type continued 
fraction 


1+ 2— 3+ 2— 5+ 2—7+ 2 
The maximum relative error of this approximation is 2-9. 
5. Multiply 2-4 by 2-°. 


6. Finally, add the hexadecimal exponent a to the characteristic of the answer. 
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Effect of an Argument Error 


«e ~ A. If the magnitude of x is large, even the round-off error of the argument 
causes a substantial relative error in the answer because A = 8 * x. 


xxxSGAMA Subprogram (GAMMA and ALGAMA) 

Algorithm 

1. If0 < x S 2-52, then compute log-gamma as log. T (x) = —log, (x). This 
computation uses the real logarithm subprogram (xxxsLoc). 

2. If 2-252 < x < 8, then compute log-gamma by taking the natural logarithm of 
the value obtained for gamma. The computation of gamma depends upon the 
range into which the argument falls. 


I (x +1) 
% 


3. If 2-952 < x < 1, thenuseI (x) = to reduce to the next case. 


4, If 1 =x S 2, then compute gamma by the following approximation: 
T (x) = do + 4% + gz? +... + 92° 

where z = x — 1.5. The coefficients were obtained by expanding the function 
f(z) = I (x) with respect to the Chebyshev polynomials for |z| = 0.5. The 
absolute error of this approximation is less than 1.5 + 2-. 

5. If2 <x < 8, thenusel (x) = (x — 1) I (x — 1) to reduce step by step to 
the preceding case. 

6. If 8 = x, then compute log-gamma by the use of Stirling’s formula: 


log, T (x) = x(log, (x) — 1) = 5 log, (x) +5 log, (Qe) + G(x). 


The modifier term G(x) is computed as 
G(x) = byx-1 + box-?. 
The absolute error of the approximation for G(x) is 1.4 * 2-28. This computa- 
tion uses the real logarithm subprogram (xxxsLoc). 
For gamma, compute I (x) = e%, where y is the value obtained for log- 
gamma. This computation uses the real exponential subprogram (xxxsExP). 


Effect of an Argument Error 


« ~ yw (x) * A for gamma, and E ~ y (x) * A for log-gamma, where y is the 
digamma function. 


If > <x < 3, then —2 < y (x) < 1. Therefore, E ~ A for log-gamma. How- 


ever, because x = 1 and x = 2 are zeros of the log-gamma function, even a small 
3 can cause a substantial « in this range. 

If the value of x is large, then y (x) ~ log, (x). Therefore, for gamma, 
e ~ 5x ° log, (x). In this case, even the round-off error of the argument con- 
tributes greatly to the relative error of the answer. For log-gamma with large 
values of x, « ~ 8. 


xxxSLOG Subprogram (ALOG and ALOG10) 


Algorithm 
1. Write x = 16? * m, where p is an integer and m is within the range, =<m<l. 
2. Define two constants, a and b, where a = base point and 2~® = a, as follows: 


1 l Bee 
lit Sm < -g thena = 7G andb = 4. 
1 1 _1 o 
lf Sm < 7, thena = ri and b = 2, 

1 


If—— =m < 1,thena =landb = 0. 


bo 











, _ m—a uote. alee ek 
3. Write z = are . Then, m = a ia ,and |z| = z° 
4. Now, x = 24? e ath and log.x = (4p — b) log.2 + logs ( a ). 


5. Finally, log. Gj = ‘) is evaluated using the Chebyshev interpolation of degree 





: 1 . ' ; 
4 in 2? over the range, OS 2? = 7° The maximum relative error of this ap- 
proximation is 2~27-8, 


6. If the common logarithm is desired, then logyox = logice * log.x. 


Effect of an Argument Error 


E ~ 8. Specifically, if 8 is the round-off error of the argument, e.g., 8 ~ 6 * 10-8, 
then E ~ 6 * 10-8. Therefore, if the argument is close to 1, the relative error can 
be very large because the value of the function is very small. 


xxxSSCN Subprogram (SIN and COS) 











Algorithm 
1. Define z = er x| and separate z into its integer part (q) and its fraction part 
T 
(7). Then z = q +1, and |x| = ( ) + ( ql +r), 


2. If the cosine is desired, add 2 to q. If the sine is desired and if x is negative, 
add 4 to q. This adjustment of q reduces the general case to the computation 
of sin (x) for x = 0 because 


cos (+x) = sin a: «) , and 
sin (—x) = sin (+ + x). 
3. Let qo = q mod 8. 


Then, for go = 0, sin (x) = sin ( i ° r) 





do = 1, sin (x) = cos Fil — ")) 
Go = 2, sin (x) = cos (7 r) 
do = 3, sin (x) = sin (F-1)) 





do = 4,sin (x) = —sin Ge r) 
go = 5, sin (x) = —cos (Fa = )) 
Jo = 6, sin (x) = —cos (+ r) 


Qo = 7, sin (x) = —sin va - )) 





. nee 
These formulas reduce each case to the computation of either sin ( a n) or 


cos (+ 1, ) where 1; is either r or (1 — r) and is within the range, 


0OfnL1. 
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4, Finally, the computation for either the sine or the cosine is performed, using 


the Chebyshev interpolation of degree 3 in 1,2. The maximum relative error 
of the sine polynomial is 2-81 and that of the cosine polynomial is 2~*+-°. 


Effect of an Argument Error 


E ~ A. As the value of x increases, A increases. Because the function value dimin- 
ishes periodically, no consistent relative error control can be maintained outside 
Tv 


the principal range, —- SxS + 9 


2 


xxxSSCNH Subprogram (SINH and COSH) 


Algorithm 


1. If |x| < 0.3465736, then compute sinh(x) as: 


sinh(x) = x + 0.16666505x? + 0.00836915x°. . 
sinh(x) 


The coefficients were obtained by expanding the function f(z) = aera 


= x, with respect to the Chebyshev polynomials over the range 
< z < 0.12011326. The relative error of this approximation is less than 
— 26.5 


Ww OX 


2. If either |x| = 0.3465736 or the cosh(x) is desired, obtain w = e'!. Then, 


w+tw} —w 
2 2 
tial subprogram (xxxsExP) is used to compute the value of w. 


af 
cosh(x) = , and sinh(x) = (sign x). = . The real exponen- 


Effect of an Argument Error 


For the hyperbolic sine, E ~ A* cosh(x) and e ~ A* coth (x). 
For the hyperbolic cosine, E ~ A sinh (x) and e ~ 8° tanh(x). 

Specifically, for the cosine, e ~ A over the entire range; for the sine, e ~ 8 for 
small values of x. 


xxxSSQRT Subprogram (SQRT) 
Algorithm 
1. Ifx = 0, then the answer is 0. 


2. Write x = 167? * m, where p is an integer and m is within the range, 
1 
cs 
356 =™ <1. 
3. Then, \/x = 16? « \/m, where p is the exponent of the answer and m is the 
mantissa of the answer. 


4. For the first approximation of \/m, take hyperbolic approximations of the form 


a+ ae where the values of a, b, and c depend upon the value of m as follows: 


c+x 
“ ts <m <l,then a = 1.80713 
b = —157727 
c = 0.954182 


These values minimize the maximum relative error (e9) over the range, 
while making an exact fit atm = 1. The exact fit at m = 1 minimizes the 
computational loss of the last hexadecimal digit for the values of m slightly 
less than 1. The relative error of this approximation is less than 2—5-*, 


be if — 2% <g ‘then a = 0.428795 
= —0.0214398 
c = 0.0548470 


These values minimize m1/8 + ¢) over this range of m where «9 denotes the 
relative error of this approximation. ¢ is less than 2—®5 * m—1/8, 
5. Multiply the result by 16? to obtain the first approximation (yo) of the answer. 
6. To obtain the final answer, the Newton-Raphson iteration 


Ynti = z( ti) 


must be applied twice to yo. For 79 <= m < 1, the final relative error is 


the final absolute error is 


‘ 1 
— 24. T. 
theoretically less than 2- for >>> 356 = =Sm<a 16” 


theoretically less than 2—?° + 16?. 


Effect of an Argument Error 


em oz 6. 


2 


xxxSTANH Subprogram (TANH) 


Algorithm 
1. If |x| = 2-%, then tanh (x) = x. 
2. If 2-1? < |x| < 0.54931, use the following fractional approximation: 
tanh (x) 1- x? + 35.1535 
- x? + 45,1842 + OO 
This approximation has a relative error less than 2—2%. The formula can be 
obtained by transforming the continued fraction 


tanh (x) _ 4g Ce i, 
x 3+ 5+ T+ w 
where w has an approximate value of 0.0307, but the true value of w is 
a 
9+ 1i+°"° 
3. If 0.54931 = x < 9.011, then use the identity tanh (x) = 1 — ey The 


computation for this case uses the rea] exponential subprogram (xxxsExP). 
4, If x => 9.011, then tanh (x) = 1. 


5. Ifx = —0.54931, then use the identity tanh (x) = —tanh (—x). 
Effect of an Argument Error 


E ~ (1 — tanh? x) A, ande ~ ee For small values of x, « ~ 8, and as the 
sinh (2x) 


value of x increases, the effect of § upon « diminishes. 


xxxSTNCT Subprogram (TAN and COTAN) 


Algorithm 


1. Divide |x| by and separate the result into the integer part (q) and the 


fraction part (r). Then, |x| = “7 (4 + r). 
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2. Obtain the reduced argument (w) as follows: 
if gq is even, then w = 1, 
if g is odd, thenw = 1 — 1. 
The range of the reduced argument is 0 = w = 1. 


tn (4 7 w) and cot |x = cot ‘s 


go = 1, tan [x = cot ( 7 ° w) and cot |x| = tan 


3. Let qo = q mod 4. 








Then, for qo = 0, tan |x| 


:) 
=) 
7 : w) and cot |x| = —tan (2 ee w) 


qo = 3, tan |x| = —tan ( 5 ° w) and cot |x| = —cot Ce w) 


7 ° w) and cot ( ji . w) are computed as the ratio of 


2 
° re )= ws Piw?) and cot G 
QO(w?) ” 
where P(w?) = ob 58037 — 12.559912w? 


QO(w?) = 270.665736 — 71.645273w? + wt 
This approximation is obtained by economizing the continued fraction 














qo = 2, tan |x] = —cot ( 














4. The values of tan 


two pene 








tan 


+o) = Ow) _Q(w?) 


we P(w?) 





tan(z) _ 1 ee ee 
Zo ae a 
in the following way: 
2 2 2 
Write: antZ ean) : eee 2 2 - 


(3 +d,)— (5+ ds)— (7+ ds) 
and determine values for d;, dz, and ds so that the right-hand expression 
gives the exact answers for z? = 0.19, 0.432, and 0.594. Then the maximum 


relative error of this formula over the range 0 S z S—~is 1.74 + 10-8. 








4 
Change the variable from z to w = ae z and rewrite the formula to ob- 
T 

tain P(w?) and QO(w?). 
5. If x < 0, then tan |x|, and cot(x) = —cot |x]. 
Effect of an Argument Error 

A 2 
eas = : oe 
en Ga? and « al oxy (Sx) for tan(x). Therefore, near the singularities 


1 ; pte 
= ( +z)» where k is an integer, no error control can be maintained. 


This is also true for cotan(x) for x near kr, where k is an integer. 


Appendix B contains accuracy and timing statistics for 
the explicitly called mathematical subprograms. These 
statistics are presented in Table 14 and are arranged 
in alphabetical order, according to the entry names. 
The following column headings are used in Table 14: 


Entry Name: This column gives the entry name that 
must be used to call the subprogram. 


Argument Range: This column gives the argument 
range used to obtain the accuracy figures. For each 
function, accuracy figures are given for one or more 
representative segments within the valid argument 
range. In each case, the figures given are the most 
meaningful to the function and range under consid- 
eration. 


The maximum relative error and standard deviation 
of the relative error are generally useful and revealing 
statistics; however, they are useless for the range of a 
function where its value becomes 0, because the slight- 
est error in the argument can cause an unpredictable 
fluctuation in the magnitude of the answer. When a 
small argument error would have this effect, the maxi- 
mum absolute error and standard deviation of the 
absolute error are given for the range. For example, 
absolute error is given for sin(x) for values of x near z. 


Sample: This column indicates the type of sample 
used for the accuracy figures. The type of sample de- 
pends upon the function and range under consid- 
eration. The statistics may be based either upon an 
exponentially (E) distributed argument sample or a 
uniformly (U) distributed argument sample. 


Accuracy Figures: This column gives accuracy fig- 
ures for one or more representative segments within 
the valid argument range. The accuracy figures sup- 
plied are based upon the assumption that the argu- 
ments are perfect (i.e., without error and, therefore, 
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having no error propagation effect upon the answers ). 
The only error in the answers are those introduced by 
the subprograms. Appendix A contains a description 
of some of the symbols used in this appendix; the 
following additional symbols are used in the presenta- 
tion of accuracy figures: 

The maximum 
relative error 


produced 
during testing. 


M (e) = Max 








The standard 
deviation (root- 
mean-square ) of 
the relative error. 


2 


1 f(a) — g(a) 
o {e) wht i fa 





The maximum 
absolute error 
produced 

during testing. 


M (E) 


Max | f(x) — g(x) | 


The standard 
deviation (root- 
mean-square) of 
the absolute error. 


E(x) — g(a)} 








o (E) = Ei 


In the formulas for the standard deviation, N repre- 
sents the total number of arguments in the sample; 
iis a subscript that varies from 1 to N. 

Accuracy for the Model 44 is based on performance 
with the FLOATING-POINT PRECISION switch in the “14” 
(vertical) position. This position selects the highest 
(56) of the four long-precision increments permitted. 

Average Speed: This column gives the timing statis- 
tics. These statistics represent the average speed in 
microseconds for the various System/360 models. Sta- 
tistics are supplied for Models 30, 40, 44, 50, 65, and 
75. Statistics for the Model 75 are based upon two-way 
interleaving. 
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Table 14. Performance Statistics 


















Accuracy Figures Average Speed ( Microseconds ) 


Entry Argument Sample 
Name Range E/U M (e) o (e) M (E) o (E) 30 | 40] 50] 65] 75 44 
(See Note 8) | (See Note 9) 


ALGAMA | 0<x <05 109x 10-6 | 3a5x10-7 | = |————S—SsY- 10800] 2800 727 553 






0B<x<3 ee Ys 734560 
iaixi0-* | 2eexu0-7 | 12100 | 3250 [1020 520 638 
8<x<16 125x10-@ | 389x10-7 | ss 7600 | 2010 162 470 356 

6<x<500 104x10-6 | 203x10-7 [| 7600 | 2010 477 362 
x<0.5,x >1.5 


l ; 

055x515 PE 846 x 10-7 | 8.62 x 10-8 4481| 1178} 361 52.3 | 229 176 
8.32x 10-7 | 120x10-7 | 361| 91.9] 52.3] 92 

ALOGIO [0.5=x<=15 Po 64K 10-7 | 478 xX 10-8 4847 | 1278 248 193 


5340} 1460] 451 | 118 70.6 325 238 


For ATAN in FORTRAN IV (E) 
3602 913 255 688 40.6 
Additional time for ATAN in 
FORTRAN IV 


x<0.5,x > 15 1.05 x 10-6 2.17 X 10-7 


: 
> 
[e.) 
pea 
iy 
q 
oo 


>| > > 
a) 8 : 
—_— 

Z19 Q 







-lSxsS+1 U 8.56 x 10-7 2.38 X 10-7 


The full range Note 7 9.75 X 10— 4.54 x 10-7 


(Mod. 44: 
The full range Note 7 9.75 X 10— 4.54 X 10-7 









ATAN 









—_+~— 


"oT 
2 2 









omic 
8 


5 
8 
2 


15783 | 4433/1329 203 1089 


14021] 3061] 638 | 150 733 65 
4 


48343 |11705 |2335 |507 | 301 3020 2807 
41737 |10144}2048 | 456 | 260 2616 2419 


CDEXP = | |x| SJ, [xl S U 4.04 x 10- 1.39 X 10-16 


tan y, y in 
CABS The full range Note 1 1.87 x 10-6 7.65 X 10-7 
7.21 X 10-7 
Ix:] S 20, |xe| S 20 U 3.63 X 10-15 | 1.29 x 10-16 


qT 
7 
[xa] = 10, |xe| — 20 U 1 79 x 10-6 
See Note 2 
CDABS Note 1 3.32 X 10-15 | 5.16 x 10—16 
16 


CDCOS [x:| = 10, |xe| =1 U 5.16 x 10—15 3.42 X 10—16 
See iJote 3 


CDLOG _} The full range Note 1 8.73 X 10-15 | 6.38 x 10-17 


CDSIN |x| S 10, [xe] S 1 U 3.72 * 10-15 3.49 x 10—16 
See Note 4 


CDSQRT | The full range Note 1 9.86 x 10-16 | 1.91 xX 10~—16 


CEXP |x:| S 170, [xe] <> U 118 x 10-6 | 2.34 x 10-7 


[xs = 170, 
6 -—7T 
£ < |xe| < 20 U 1.06 x 10: 2.51 X 10 
CLOG The full range Note 1 2.00 x 10-6 1.56 X 10-7 
except (1 + Oi) 


53362 }11940}2393 | 542 | 317 3105 2860 


6 
6 


Q 
n 4 
bo 


48250 {11668} 2322} 503 j 298 3018 2808 


27951 | 5906/1282; 301 | 181 1543 1392 


0 
4 
13731} 3888} 1166} 291 | 177 958 73 
4 


13899} 3930] 1180] 294 | 178 
15504| 4246] 1261] 338 | 216 1014. 777 
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Table 14. Performance Statistics (Continued ) 


Accuracy Figures Average Speed ( Microseconds ) 
E/U 


M (e) o (e) M (E) o (E) 30 40 | 50} 65 | 75 44 
(SeeNote8) |(See Note 9) 
303 


























1.35 x 10-7 


U 
U 
U 
U 
3.80 X 10-4 7.70 X 10-4 4610 ; 1220} 351) 89.3 | 55.5 
U 
See Note 5 
U 1.13 X 10-5 6.03 xX 10-7 4580 | 1210 | 348 | 88. 
See Note 5 
U 1.67 x 10-5 6.67 X 10-7 
See Note 5 
CSIN |x:| S 10, |xe| S1 U 1.97 x 10-6 7.09 x 10-7 
See Note 6 
2 





10 < |x| S 100 


Ot gt 
OL OL 
i) oO 


o> 
— 
i) 
S 
Q 
— 
S 
oe] 
— 
i?) 
> 
i) 


15690 | 4397 |1316 | 33 


28701 805| 219 676 493 


For DATAN in FORTRAN IV (E) 
19056 4000 715 153 83.6 
Additional time for DATAN in 

FORTRAN IV 
104 


8 


5.93 X 10-17 13133 | 3146} 605} 132 


1.01 x 10-15 13133 | 3146 


15300 | 3590 | 67 
15900 } 3640] 71 
15800 | 3690 | 70 


15800 | 3690} 70 


CSQRT The full range Note 1 1.61 x 10-6 4.58 X 10-7 
DARCOS 2.72 x 10-16 | 9.35 x 10-17 
DARSIN |-1Sx<S +1 


DATAN The full range Note 7 2.08 X 10-16 | 6.64 x 10-17 
(Mod, 44: 
tan y, y in 
eis Bete: 
2, 
v¢ 


2.40 X 10-16 6.00 x 10—17 










DATAN2 | The full range Note 7 2.08 x 10-16 | 6.64 x 10-17 


-—-10sx<0 —16 


: 


DCOTAN 3.46 x 10-16 | 8.38 x 10-1 
Hig ee 1.72 x 10-18 | 5,00 x 10-15 
4 = See Note 5 
7 U 5.33 X 10-18 | 1.09 x 10—14 
Zsmls | OY | sees | 


10< |x| = 100 U 8.61 X 10-18 4.61 X 10-14 
See Note 5 
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Table 14. Performance Statistics (Continued ) 


Accuracy Figures Average Speed ( Microseconds ) 
E/U 


M (e) o (e) M (E) o (E) 30 40 | 50] 65 | 75 44 
(See Note 8) | (See Note 9) 


perF [pjsisi |u| itoxio-| amxi-m | | 38400] 010 1208 1083 












Argument 
Range 





q 










q 










q 


bo 


45200 }10900 |2090 2871 2628 


12145 2907| 607) 648 
12145| 2907| 607 716 644 


12145 | 2907} 607 : 715 


| 
bo 
o 


DERFC [5 ser00 | 6920 [1700] 360 ax8 | 3055 s1a7 
18600 12131088 
1.317 <x 2.04 24100 | 6000 [1130] 244 [139 | 1623 1454 
2.04<x<4 45000 [10800 [2060] 444 [258 | 2854 2612 

i: 

; 

: 


PSI 
ee ee es 
cw] wt Ot 
Co} Cj = 
q |S 
1; 
Su] ou] or] Ot] G@ 
Out out Ot 


— 
oy) 
ee) 
fon) 
y 
eo 


DGANMA 


F-3i98 —3204 | 


ie) 
ho 
es 
i" 


2498 2294 
4] 2503 2298 
3056 2783 
2983 2709 
3228 2947 
1712 1572 
1725 1585 


6.42 X 10-15 } 2.01 x 10-15 40600 | 9560 |194 


iw) 
oe | a 
< (|S 
bo 


ae 
nse 
ol1o 


62x 10-14 | 2.96 x 10-14 40600 | 9560 |19 
4.11 x 10-36 | 1.60 x 10-16 | sf 46900 |11300 216 


16x <0 57 
DLGAMA |0<x<05 
05<x<3 Pf 2.86 X 10-16 | 1.16 x 10-18 44900 {11100 
j 2.38 X 10-15 | 3.99 x 10-16 45900 {12200 
6 
| 


| 
— 
o> | Go 
S/S 
iar 
3 
S| ia 
bo | po 
ona es) 
rsa RS 


be 
bdo 
i<o) 
bo 


dO | Ot 
co 
oo a) 
fm 
on) 
— 


an ae 1.85 X 10-16 16044] 3769] 734 


92.0 
4.08 X 10-16 | 4.85 x 10-17 9.10 X 10-17 2.17 X 10-17 13145 | 3148 } 609] 133 | 755° 762 714 
a 1.64 X 10—16 6.35 X 10-17 13145} 3148 } 609] 133 75.5 763 715 


2.69 X 10-15 1.03 X 10-15 13145 | 3148 | 609] 133 | 75.5 1 761 714 
5.29 X 10-17 86501 2170] 408| 88.5 505 457 
19 


3.59 x 10-16 | 8.73 x 10-17 18400 | 4280 207 1049 925 


No 
ico} 
jon) 
—_ 
o> 
—_ 






pt 
o>) 
— 
oo 
[o>] 
on) 
co 
it) 
co 
oO 
OL 
Ot 


(.) 
jor) 
ron 


929 855 


(oy) 


CO 
ll Bleed lad 
o|l=| = 

co 

bo 


< 
DLOG 052x215 
x<05,x > 1.5 
DLOGIO |05Sx15 
a 


oo 


bo 


x<0.5,x > 1.5 
DSIN kj <= 
2 
[ty </l S10 
10 < |x|S 100 


DSINH | |x| S$ 0.34657 
0.34657 < |x| $5 
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Table 14. Performance Statistics (Continued ) 


Entry 
Name 
See Note 5 


DSQRT The full range E 1.08 X 10-16 | 2.17 xX 10-17 
1.57 X 10-18 


DTAN U 5.25 x 10-16 | 9.96 x 10-27 
pea Penal U 
4 <ixis 5 
= < |x| <10 
2 _ See Note 5 



































Accuracy Figures Average Speed (Microseconds ) 


M (E) o (E) 30 | 40 | 50| 65 | 75 4A 
(See Note 8) | (See Note 9) 


a ri fo 


12299 | 2850 | 477/106 55. 
16078 | 3778 | 833) 19% 110 


4183 | 1250 
4500 | 1230 | 397 
10000 | 2800 | 845/213 [127 
7040 | 1960 
4250 | 1150 | 9s e 
9950 8351210 
10100 | 2840 | 859/216 
5840 484|123 o{ 4 
5620 489|12 73.5 | 4 


Argument Sample 


Range M (e) o (e) 










1.67 X 10-12 | 3.69 x 10-14 


q 


ty 


4.51 xX 10-15 


1lo< |x| = 100 U 3.79 X 10-12 | 9.50 x 10-14 

See Note 5 
DTANH |x| S 0.54931 2.00 x 10-16 4.45 x 10-17 
0.54931 < |x| $5 1.99 x 10—16 9.54 x 10-17 


OU 


to 


8 

6 

9 

4.65 xX 10-7 1.28 x 10-7 3 
4.69 xX 10-7 1.17 x 10-7 
9.26 X 10-7 1.43 x 10-7 

9.02 x 10-8 | 3.42 x 10-8 8 
6.07 x 10-8 3.42 x 10-8 
9.10 x 10-7 | 2.97 x 10-7 

3.90 x 10-6 | 5.65 x 10-7 2 

1.02 x 10-6 | 213 x 10-7 3 

1.20 x 10-6 | 3.60 x 10-7 

152x 10-5 | 845 x 10-6 

9.86 x 10-7 | 3.45 x 10-7 

1.00 X 10-7 | 3.74 x 10-8 


|x| = 170 


ERF |x| S 1.317 
1.317 < |x| S 2.04 
2.04 < |x| S 3.9192 


ERFC -38<x<0 
0=x£1.317 
1.317 <x = 2.04 
204<x<4 
42x2133 

GAMMA |0<x<l 
lsx22 
9<x=4 
4<x<8 
8=x<16 
16=x<57 

SIN 


Ix| 


We) 
n= 
Ot 
oe) 
aN 


we 
“l 
jor) 


ie) 
— 
wt 
ie) 
ee) 


~] 
aS 
Oo 


jon] 


3 
5 


jor) 
an) 
~l 
_ 
1) 
(eS) 
co 
—_ 
fon) 


(ew) 
Pata} 
tN 
Ko) 
a, we 
o!]1o 
Ot 
nS 


on) 
~l 
ow 
1 


~l 
cS 
i<e) 


9 
6 
7 
1 
0 
2 
2 
8 
] 
3 


~1 
aN 





eX) 


9.29 x 10-7 | 3.63 x 10-7 6330 |1700 | 546/137 | 81.1 
2.95 x 10-6 | 8.14 x 10-7 7740 |2070 | 659|166 | 96.2 


2.29 X 10-5 
4.36 x 10-5 


7.67 X 10-6 
1.45 x 10-5 


_ 
oO 
No 
oO 
bo 
fo») 
(ev) 


[ee] 
_ 
: 


12000 | 3320 


12000 | 3320 110201263 8 


3876 } 1036 } 298} 74.2 | 44.3 19 
3980 | 1064 201 163 


| 210 142 90 


oO 
ro 
ine) 


3 846 
8 641 
9 874 
1 192 
8 189 
189 
2 220 
525 
1 367 
4 195 
9 216 
521 
532 
3 306 
301 
330 
402 
618 
616 
4 155 


1.59 x 10-6 2.02 X 10-7 1.31 X 10-7 5.55 X 10-8 


5.53 X 10-8 


bo] 3 

| MA 
* jtols 
A 

= 


1.41 X 10-7 









10 < |x| = 100 
-5Sx2+5 


The full range 


1.46 x 10-7 5.61 xX 10-8 3989 | 1064 
5890 | 1740 


2965 


1.20 x 10-6 
8.70 X 10-7 


3.20 X 10-7 


E 1.68 X 10-7 





ad 
fon) 


Table 14. Performance Statistics (Continued ) 









Accuracy Figures Average Speed ( Microseconds ) 


7 
E/U ete) M (E) o (E) 30 | 40 | 501 65 | 75 44 
(See Note 8) | (See Note 9) 
> <he <_7 6.58 x 10-5 1.67 x 10-6 4500 | 1184} 338 | 85.3 | 52.8 
|S5 See Note 5 






Argument 


4,92 x 10-5 1.28 x 10-6 

7 < 

eos 3.35 X 10-5 1.02 x 10-6 4460 |1170; 335 | 84.1 | 52. 
See Note 5 


assert a ca 
eaRE | U | «5.74 x 10-7 7.53 x 10-8 ior Paes: 5052 551 (142 | 86.3 


Notes to Table 14 

These notes are associated with Table 14 and contain 
more detailed information about samples and relative 
errors for certain functions. 


Note 1: The distribution of sample arguments upon 
which these statistics are based is exponential radially 
and is uniform around the origin. 


Note 2: The maximum relative error cited for the 
ccos function is based upon a set of 2000 random argu- 
ments within the range. In the immediate proximity of 


the points (» + +) a + 0i (where n= 0, +1, +2, 


,) the relative error can be quite high, although 
the absolute error is small. 


Note 3: The maximum relative error cited for the 
cpcos function is based upon a set of 1500 random 
arguments within the range. In the immediate prox- 


2 
+1, +2,...,) the relative error can be quite high 


imity of the points (» + +) a + Oi (where n = 0, 


although the absolute error is small. 





Note 4: The maximum relative error cited for the 
cpsin function is based upon a set of 1500 random 
arguments within the range. In the immediate prox- 
imity of the points nx + Oi (where n = +1, +2, 

,) the relative error can be quite high although 
the absolute error is small. 


Note 5: The figures cited as the maximum relative 
errors are those encountered in a sample of 2500 ran- 
dom arguments within the respective ranges. See the 
appropriate section in Appendix A for a description 
of the behavior of errors when the argument is near a 
singularity or a zero of the function. 

Note 6: The maximum relative error cited for the 
csin function is based upon a set of 2000 randon argu- 
ments within the range. In the immediate proximity 
of the points nx + 0i (wheren = +1, +2,...,) the 
relative error can be quite high although the absolute 
error is small. 


Note 7: The sample arguments were tangents of 


numbers uniformly distributed between — a and + De 


Note 8: The statistics for the Model 75 are based 
upon two-way interleaving. 


Note 9: The second column of speeds for the Model 
44 applies to that machine with high-speed registers. 


227 =6184 
89 63 
446 294 


Appendix C. 


Appendix C contains descriptions of the procedures 
followed when the execution of a load module is dis- 
continued. Execution may be discontinued due to one 
of two reasons: an interruption or an error. After an 
interruption is processed, execution of this load module 
or phase continues; after an error is processed, execu- 
tion of this load module or phase is terminated. The 
following text explains the procedure used to handle 
each case. 


Interruption Procedures 


A program interruption is a computer-originated 
break in the flow of processing. (For a full description 
of program interrupts, see the publication IBM Sys- 
tem/360 Principles of Operation, Form A22-6621.) The 
FORTRAN library processes those interrupts that are 
described below; all others are handled directly by 
the system Supervisor and will cause job termination. 


The following services are provided by the rorTRAN 
support for interrupt occurrence: 


1. When an interrupt occurs, indicators are set 
to record exponent overflow, underflow, fixed-point, 
floating-point or decimal divide exceptions. These 
indicators can be interrogated dynamically by the sub- 
programs described in “Service Subprograms.” 


2. A message is printed on the object error unit 
when each interrupt occurs. The psw printed in the 
message indicates the cause of each interrupt. 


3. Result registers are changed when exponent over- 
flow or exponent underflow (codes C and D, below) 
occur. (For a description of the format of floating- 
point numbers, see the publication IBM System/360 
Principles of Operation, Form A22-6621.) Result regis- 
ters are also set when a floating-point instruction is 
referenced by an assembler language execute (Ex) 
instruction. 


IHC2101 PROGRAM inTeRRuPT ( . Jou PSW IS xxxxxxx 


4 
5 
6 
7 
9 
Cc 
D 
F 





Figure 1. Format of Program Interrupt Message, Operating System 


Interruption and Error Procedures 


4, Condition codes set by floating-point addition or 
substraction instructions are altered for exponent 
underflow (code D). 


5. After the foregoing services are performed, execu- 
tion of the program continues from the instruction 
following the one that caused the interrupt. 

The program interrupt message contains the old pro- 
gram status word (psw), which indicates the cause 
of the interrupt. Figure 1 shows the format of the 
message as it is issued by the operating system. Figure 
2 shows the format as issued by the Model 44 system, 


| and Figure 2.1 shows the pos format. If the letter A 


appears in parentheses in the program interrupt mes- 
sage, boundary adjustment has taken place. The letter 
P in the message indicates that the interruption was 
precise. This will always be the case for nonspecifica- 
tion interrupt messages in FORTRAN except when using 
machines with special hardware on which imprecise 
interruptions may occur. The eighth character in the 
psw represents the code number associated with the 
type of interruption. These interruptions are described 
in the following paragraphs. (For more information 
on the psw, see the publication IBM System/360 Prin- 
ciples of Operation, Form A22-6621. ) 


Specification Exception (Code 6): The specification 
exception (code 6) is recognized when a data address 
does not specify an integral boundary for that unit of 
information. A specification error would occur, for ex- 
ample, during the execution of the following program 
segment: 


DOUBLE-PRECISION D,E 
COMMON A,B,C 
EQUIVALENCE (B,D) 

D = 3,0D02 


Fixed-Point Divide Exception (Code 9): The fixed- 
point divide exception (code 9) is recognized when 
division of a fixed-point number by zero is attempted. 


XXXXXXXX 
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A fixed-point divide exception would occur during exe- 
cution of the following statements: 


J=0 
I=7 
K=I/J 


Exponent-Overflow Exception (Code C): The ex- 
ponent-overflow exception (code C) is recognized 
when the absolute value of the result of a floating- 
point addition, subtraction, multiplication, or division 
is greater than or equal to 16° (approximately 7.2 x 
1075). For example, an exponent overflow would occur 
during execution of the statement: 


A= 1.0E + 75 + 7.2E + 75 


When the interrupt occurs, the result register con- 
tains a floating-point number whose fraction and sign 
is correct. However, the number is not usable for 
further computation since its characteristic field no 
longer reflects the true exponent. The content of the 
result register as it existed when the interrupt occurred 
is printed following the program interrupt message 
with the format. 


REGISTER CONTAINED hhhhhhhhhhbhhhhhh 
where: 


hhhhhhhhhhhhhhhh is the floating-point 


number in hexadecimal 
notation. 


Exponent overflow causes “exponent wraparound” 
—i.e., the characteristic field represents an exponent 
that is 128 smaller than the correct one. Treating bits 1 
through 7 (the exponent characteristic field) of the 
floating-point number as a binary integer, the true ex- 
ponent (TE) may be computed, as follows: 


TE = (Bits 1 through 7) + 128 — 64 


Before program execution continues, the FORTRAN 
library sets the result register to the largest possible 
floating-point number that can be represented in 
short precision [16® * (1-16—®)] or in long precision 
[166° * (1-16—14)], but the sign of the result is not 
changed. The condition code is not altered. 

Exponent-Underflow Exception (Code D): The ex- 
ponent-underflow exception (code D) is recognized 
when the absolute value of the result of a floating-point 
addition, subtraction, multiplication, or division, is less 
than 16—® (approximately 5.4 x 10-7°) but not equal 
to 0. An exponent-underflow exception would occur 
during execution of the statement: 

A=-10E-50 * 10E—50 

Although exponent underflows can be masked, 
FORTRAN jobs are executed without the mask so that 
the library will handle such interrupts. 

When the interrupt occurs, the result register con- 
tains a floating-point number whose fraction and sign 
is correct. However, the number is not usable for fur- 


48 


ther computation since its characteristic field no longer 
reflects the true exponent. The content of the result 
register as it existed when the interrupt occurred is 
printed following the program interrupt message with 
the format: 


REGISTER CONTAINED hhhbhhhhhhhbhhhh 


where: 


hhhhhhhhhhhhbhhh is the floating-point 
number in hexadecimal 
notation. 


Exponent underflow causes “exponent wraparound” 
—i.e., the characteristic field represents an exponent 
that is 128 larger than the correct one. Treating bits 1 
to 7 (the exponent characteristic field) of the floating- 
point number as a binary integer, the true exponent 
(TE) may be computed, as follows: 


TE = (Bits 1 through 7) — 128 — 64 


Before program execution continues, the FORTRAN 
library sets the result register to a true zero of correct 
precision. If the interrupt resulted from a floating-point 
addition or subtraction operation, the condition code 
is set to zero to reflect the setting of the result register. 
Note: The Sysem/360 Operating System FORTRAN pro- 
grammer who wishes to take advantage of the “ex- 
ponent wraparound” feature and handle the interrupt 
in his own program must call an assembly language 
subroutine to issue a SPIE macro-instruction, which will 
override the FORTRAN interruption routine. 

Floating-Point Divide Exception (Code F): The 
floating-point divide exception (code F) is recognized 
when division of a floating-point number by zero is 
attempted. A floating-point divide exception would 
occur during execution of the following statements: 


B = 0.0 
A= 1.0 
C = A/B 


System/360 Operating System 


The interrupt message is written in the system out- 
put data set. 

A specification exception program interrupt message 
(code 6) is issued only if the BoUNDRY= ALIGN option 
was specified in the FoRTLIB macro-instruction during 
system generation and a boundary alignment error 
occurs. Then the boundary alignment routine is in- 
voked to correct the boundary misalignment. If an 
instruction that has been processed for boundary mis- 
alignment also contains a protection, addressing, or 
data error, the interrupt message will be reissued with 
the appropriate code (4, 5, or 7). (In these cases, the 
letter A appears in parentheses in the program inter- 


OA2101 


PROGRAM INTERRUPT (A) OLD PSW JS xxxxxxx 


Figure 2. Format of Program Interrupt Message, Model 44 System 


rupt message.) Then the job will terminate because 
both a specification error and a protection, addressing, 
or data error have been detected. The completion code 
in the dump will specify that the job terminated be- 
cause of the specification error. 

The number of warning messages printed is limited 
to ten. After ten boundary alignment adjustments have 
been made, the message is suppressed, but boundary 
alignment violations continue to be corrected. 

Protection Exception (Code 4): The protection ex- 
ception (code 4) is recognized when the key of an 
operand in storage does not match the protection key 
in the psw. A message is issued only if a specification 
exception (code 6) has already been recognized in the 
same instruction. Otherwise, the job terminates abnor- 
mally without a message. 

Addressing Exception (Code 5): The addressing ex- 
ception (code 5) is recognized when the address of 
the data is outside of available storage for the particu- 
lar installation. A message is issued only if a specifica- 
tion exception (code 6) has already been recognized 
in the same instruction. Otherwise, the job terminates 
abnormally without a message. 

Data Exception (Code 7): The data exception (code 
7) is recognized when the sign or digit codes for a 
CONVERT TO BINARY instruction are incorrect. A message 
is issued only if a specification exception (code 6) 
has already been recognized in the same instruction. 
Otherwise, the job terminates abnormally without a 
message. 


Model 44 Programming System 


The interrupt message is written in sysopT. 
The program interrupt message (with code 6, as 
described in “Specification Exception”) that results 


ILF2251 PROGRAM INTERRUPT — OLD PSW IS xxxxxxx 


4 





XXXXXXXX 


when the boundary specification convention is vio- 
lated will contain the “(A)” only if the «rrx option 
has been turned on (sETA 1) in BoAvopT. With that 
option on (as it is in the distributed version of the sys- 
tem), a routine that adjusts for the misalignment is 
executed each time such a violation occurs, and proc- 
essing continues. 

With the «r1x option on, the number of program 
interrupt messages put out as the result of boundary 
violations is limited to a message for each of the first 
n violations per execution, where n is the operand of 
the sera instruction for the &PRNTMEs option. In the 
system as distributed, this operand is equal to 0. Only 
the message, and not the alignment correction, is 
inhibited. 


Disk Operating System 
The interrupt message is written on SYSLST. 

The ILF225I message is issued for arithmetic excep- 
tions, such as fixed point divide. The REGISTER 
CONTAINED portion is issued only for exponent over- 
flow and underflow conditions. Processing resumes 
after the error condition is listed. 

The ILF226I message appears for other forms 
of program interrupt. The A character appears only 
for boundary alignment errors. In such cases, it is for 
information only since an attempt is made to correct 
the error and resume processing. If, however, an in- 
struction containing a boundary alignment violation 
also contains a protection, addressing, or data error, 
the message is reissued with the appropriate code, 4, 
5, or 7. The job then terminates. 

The interruption codes and their explanations are 
identical to those for the Operating System, as dis- 
cussed previously in this section. 


xXXXxxxxX [REGISTER CONTAINED xxxxxxxx] 


ILF2261 PROGRAM INTERRUPT [A] — OLD PSW IS xxxxxxx iH XXXXXXXX 
7 





e@ Figure 2.1. Format of DOS Program Interrupt Messages 
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Error Procedures 


During execution, the mathematical subprograms as- 
sume that the argument(s) is the correct type. No 
checking is done for erroneous arguments (i.e., the 
wrong type, invalid characters, the wrong length, etc. ); 
therefore, a computation performed with an erroneous 
argument has an unpredictable result. However, the 
nature of some mathematical functions requires that 
the input be within a certain range. For example, the 
square root of a negative number is not permitted. 
If the argument is not within the valid range given in 
Tables 2 through 6, an error message is written on the 
object error unit data set defined by the installation 
during system generation. The execution of this load 
module or phase is terminated and control is returned 
to the operating system. 


The error message that is issued has the format: 
IHCyyyI [Message text] 
TRACEBACK FOLLOWS... 
or 
OAyyylI 
or 
ILFyyylI 
TRACEBACK FOLLOWS... 


where yyy is a numeric code that identifies the error 
detected. 


The first message is issued by the Operating System, 
the second by the Model 44 Programming System, and 
the last by the Disk Operating System. Traceback is 
a diagnostic tool for FORTRAN under the Operating Sys- 
tem and the Disk Operating System. It is a list of rou- 
tines in the direct line of call to the routine in which 
the error occurred. It is described in the publications 
IBM System/360 Operating System FORTRAN IV 
Programmer’s Guide: FORTRAN IV (G), Form C28- 
6639, or FORTRAN IV (H), Form C28-6602, and IBM 
System/360 Disk Operating System: FORTRAN IV 
Programmer's Guide, Form C28-6397. The following 
text lists the error messages in numeric order, explains 
the error, and indicates what action the system takes. 
In these explanations, x represents the argument sup- 
plied by the programmer and * represents the largest 
possible value that can be represented in floating-point 
notation. 


1HC2161 SLITE-SLITET X IS AN ILLEGAL VALUE 

OA216I 

ILF224!1 
Explanation: In the xxx¥sLit subprogram, a value of 
i that is not 0, 1, 2, 3, or 4 is an error; for xxxsLITE, a 
value for i that is not 1, 2, 3, or 4 is an error. 
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System Action Without Extended Error Facility: Ex- 
ecution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the error 
occurrence counts kept in the Option Table. Stand- 
ard corrective action is taken before execution con- 
tinues and consists of no action for sLITE, and 
returning an OFF indication (ie., J = 2) for sLITET. 


1HC2411 FIXPI INTEGER BASE = 0, INTEGER EXPONENT = 
X, LEO 
OA2411 
| 1LF2411 
Explanation: In the xxxrixp1 subprogram, a base 
number of zero and an exponent = 0 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


1HC242I FRXPI REAL*4 BASE = 0.0, INTEGER 
EXPONENT = X, LE 0 
OA242I 
| 1LF2421 

Explanation: In the xxxFrrxpi1 subprogram, a base 
number of zero and an exponent < 0 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


1HC2431 FDXPI REAL*8 BASE = 0.0, INTEGER 
EXPONENT = X, LE O 
OA2431 
| 1F2431 

Explanation: In the xxxrpxp1 subprogram, a base 
number of zero and an exponent S 0 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


1HC244I FRXPR REAL*4 BASE = 0.0, REAL*4 
EXPONENT = X.X, LE 0 
OA244I 
| MLF2441 
Explanation: In the xxxrrxpr subprogram, a base 
n ent <= 0 is an error. 

pumiber OF aoe ane ee == 0 is ‘Te JHC2511 SQRT NEGATIVE ARGUMENT = X 
System Action Without Extended Error Facility: 6 ,o511 
Execution of this load module or phase is terminated. | F251! 


System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


Explanation: In the xxxssgrt subprogram, a value 
of x < Ois an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 


to 0. cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to |x|*. 


1HC2451 FDXPD REAL*8 BASE = 0.0 REAL*8 EXPONENT = 
X.X, LEO 
OA2451 
| F245! 
Explanation: In the xxxrpxpp subprogram, a base —_jyyje2591 EXP ARG = X.X, GT 174.673 
number of zero and an exponent <= 0 is an error. OA252! 
System Action Without Extended Error Facility: | "F252I 


Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


Explanation: In the xxxsExp subprogram, a value of 
x > 174.673 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


| 1HC246!1 FCXPI] COMPLEX*8 BASE = 0.0 + 0.01, to * 
INTEGER EXPONENT = X, LE 0 
OA246I 
| ILF2461 1HC2531 ALOG-ALOG10 ARG = X.X, LE = 0 
Explanation: In the xxxrcxp1 subprogram, a base Wirt 


number of zero and an exponent = 0 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


Explanation: In the xxxsLoc subprogram, a value of 
x < 0 is an error. Because this subprogram is also 
called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real base to a power. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 


Standard corrective action is taken before execution 
continues and consists of setting the penctions) value 
to * if x = 0, and to log |x| or logio |xp (depending on 
the subprogram called) if x is less than 0. 


| 1HC2471 FCDXI COMPLEX*16 BASE = 0.0 + 0.01, 
INTEGER EXPONENT = X, LE 0 
0A2471 
| 1LF2471 


Explanation: In the xxxFcpx1 subprogram, a base 


number of zero and an exponent = 0 is an error. IHC254I SIN-COS /ARG/ = /X.X (HEX = X)/, GE PI#2**18 
System Action Without Extended Error Facility: |OA254I 

Execution of this load module or phase is terminated. | /LF254I 

System Action With Extended Error Message Fa- Explanation: In the xxxsscn subprogram, a value of 
cility: Execution continues contingent upon the |x|= 218 « 7 is an error. 
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System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


to \/2/2. 


1HC255!1 ATAN2 ARGUMENTS = 0.0 
OA2551 
| F255! 

Explanation: In the xxxsatN2 subprogram when 
entry name ATAN? is used, a value of x; = x, = 0 is 
an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


IHC2561 SINH-COSH /ARG/ = X.X/, GE 174.673 
OA256I 
| ILF2561 

Explanation: In the xxxsscnu subprogram, a value 
of |x| = 174.673 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *. 

1HC2571 ARSIN-ARCOS /ARG/ = /X.X/ GT 1 

OA257! 

| F257! 

Explanation: In the xxxsascn subprogram, a value 
of |x| > 1 is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


IHC2581 TAN-COTAN /ARG/ = /X.X (HEX = X)/, GE 
PIe2**7 8 
OA258!I 
| F258I 
Explanation: In the xxxsrncr subprogram, a value of 
|x| = 218 + + is an error. 
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System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to l. 


IHC2591 TAN-COTAN /ARG/ = /X.X (HEX = X)/, 


APPROACHES SINGULARITY 


OA2591 
| 1LF2591 


Explanation: In the xxxstNcr subprogram, a value of 


y = 3 
x too close to one of the singularities ( + oe ot on rer 


for the tangent; +7, +2z2,... for the cotangent) is 
an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution. 
continues and consists of setting the functional value 
to *. 


1HC2611 DSQRT NEGATIVE ARGUMENT = X.X 
OA261! 
| ILF2611 


Explanation: In the xxxsgrt subprogram, a value 
of x < Ois an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to |x|*, 


IHC262I DEXP ARG = X.X, GT 174.673 
OA2621 
| ILF2621 


Explanation: In the xxxLExp subprogram, a value of 
x > 174.673 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *. ; 


1HC2631 DLOG-DLOG10 ARG = X.X, LE ZERO 
OA263! 
| ILF2631 

Explanation: In the xxxtLoc subprogram, a value of 
x <= 0 is an error. Because this subprogram is also 
called by an exponentiation subprogram, this mes- 
sage also indicates that an attempt has been made 
to raise a negative real number to a power. 


error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *, 


1HC2671 DARSIN-DARCOS /ARG/ = /X.X/, GT 1 
OA2671 
| ILF2671 


System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to — * if x = 0 and to log |x| or logio |x| (depending 
on the subprogram called ) if x is less than 0. 


Explanation: In the xxxLascn subprogram, a value 
of |x| > 1 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


1HC264I DSIN-DCOS /ARG/ = /X.X (HEX = X)/, 


* 
Sheree | 1HC2681 DTAN-DCOTAN /ARG/ = /X.X (HEX = X)/ GE 


| oe PI* (2%*50) 
OA268! 
Explanation: In the xxxiscn subprogram, a value of | F268! 


|x| = 25° » = is an error. . 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


to 2/2. 


Explanation: In the xxxLtNcr subprogram, a value of 
|x| 2 25° « wis an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 1. 


1HC2651 DATAN2 ARGUMENTS = 0.0 


rere IHC2691 DTAN-DCOTAN /ARG/ = /X.X (HEX = X)/, 
APPROACHES SINGULARITY 

Explanation: In the xxxLarne subprogram when  0A2691 

entry name DATAN? is used, a value of x, = x2 = 0 | 1F2691 

is an error. Explanation: In the xxxLtNcr subprogram, a value of 

S ° . Fe : . Ae 3 
ystem Action Without Extended Error Facility x too close to one of the singularities ( + ae + alt al 

Execution of this load module or phase is terminated. 2 270° 

System Action With Extended Error Message Fa- for the tangent; +7, +27,... for the cotangent) is 


cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0. 


an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 


IHC266I DSINH-DCOSH /ARG/ = /X.X/, GE 174. 673 
OA266! 
| ILF266! 


Explanation: In the xxx_scnu subprogram, a value 

of |x] = 174.673 is an error. IHC2711 CEXP REAL ARG = X.X (HEX =X), GT 174.673 
System Action Without Extended Error Facility: | OA2711 

Execution of this load module or phase is terminated. | /LF271! 

System Action With Extended Error Message Fa- Explanation: In the xxxcsexp subprogram, a value 
cility: Execution continues contingent upon the of x; > 174.673 is an error. 


Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *, 
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System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to * (cos x — isin x) where x is the imaginary part of 
the argument. 


1HC2751 CSIN-CCOS/IMAG ARG/ = X.X (HEX = X)/, GT 


174.673 


OA275I 
| ILF2751 


Explanation: In the xxxcsscn subprogram, a value 
of |x2| > 174.673 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 


error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 


| 12721 CEXP IMAG ARG = X (HEX = X), ABS VALUE GE continues and consists of setting the functional value 


*Q%% 

0A2721 eran as follows: If imaginary part > 174.673 (x is real 
* 

Fea part of argument): for sine, result = = (sin x + 
Explanation: In the xxxcsExp subprogram, a value 2 


of |x2| => 218 « + is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to 0 + Oi. 


* 
icos x); for cosine, result = z (cos x — isin x). If 


imaginary part < — 174.673 (x is real part of argu- 


* 
ment): for sine, result = 3 (sin x — icos x); for 


* 
cosine, result = 3 (cos x + isIN x). 


IHC2811 CDEXP REAL ARG = X.X (HEX = X) GT 174.673 
OA2811 
| ILF2811 


Explanation: In the xxxcLexp subprogram, a value 
of x, > 174.673 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to * (cos x + isin x) where x is the imaginary part 
of the argument. 


| 1HC2731 CLOG ARGUMENT = 0.0 + 0.01 
OA2731 

| ILF2731 
Explanation: In the xxxcstoc subprogram, a value 
of x, = xX. = Ois an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to — * + 04. 


IHC2821 CDEXP IMAG ARG = X.X (HEX = X) ABS VALUE 
GE PI*2**50 
OA282I 
| ILF2821 
Explanation: In the xxxcLexp subprogram, a value 
of |x| = 25° + a is an error. 
System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
Execution of this load module or phase is terminated. continues and consists of setting the functional value 
System Action With Extended Error Message Fa- to 0 + 0i. 
cility: Execution continues contingent upon the | IHC283I CDLOG ARGUMENT = 0.D0 + 0.D0! 
error occurrence counts kept in the Option Table. 0A283/ 
Standard corrective action is taken before execution | MF283I 


1HC274I CSIN-CCOS/REAL ARG/ = /X.X (HEX = X), GE 
PIE2**18 
OA274I 
| 1LF2741 
Explanation: In the xxxcsscn subprogram, a value 
of [x1 = 238 ¢ + is an error. 
System Action Without Extended Error Facility: 


continues and consists of setting the functional value 
to 0 + Oi. 


Explanation: In the xxxcLLoc subprogram, a value 
of x; = x. = 0 is an error. 


System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to — * + 0i. 


System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *. 


IHC2911 ALGAMA ARG = X.X (HEX = X), LEO OR 
GE 4.2937*10**73 
OA2911 
| fLF2911 


Explanation: In the xxxscAMa subprogram for the 


IHC284I CDSIN-CDCOS /REAL ARG/ = /X.X (HEX = X)/, 
GE PI*2**50 
OA284I 
| 1F284) 


Explanation: In the xxxcLscn subprogram, a value 
of [x1] => 250 + is an error. 
System Action Without Extended Error Facility: 


log-gamma function, a value of x S 0 or x = 4.2937 
¢ 10% is an error. 
System Action Without Extended Error Facility: 


Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


to0 + 0i. to *. 
IHC285I CDSIN-CDCOS /IMAG ARG/ = /X.X (HEX = X)/, IHC3001 DGAMMA ARG = X.X (HEX = X), LE 2** —252 OR 
GT 174.673 GE 57.5744 
OA285! OA3001 
| ILF2851 | ILF300! 


Explanation: In the xxx.camMa subprogram for the 
gamma function, a value of x = 2-8? or x = 57.5744 
is an error. ; 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 


Explanation: In the xxxcLscn subprogram, a value 
of |x2| > 174.673 is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 
Standard corrective action is taken before execution 
continues and consists of setting the functional value 
as follows: If imaginary part > 174.673 (x is real 


* to *, 
part of argument): for sine, result = 3 (sIn x + 
: : * : 1HC3011 DLGAMA ARG = X.X (HEX = X), LEO OR 
icos x); for cosine, result = 5 (cos x — isin x). If GE 4.2937*10**73 : 
imaginary part < — 174.673 (x is real part of argu- eat 
* 
ment): for sine, result = 5 (sin x — icos x); for Explanation: In the xxx~camMa subprogram for the 


log-gamma function, a value of x = 0 or x = 4.2937 
* 10% is an error. 

System Action Without Extended Error Facility: 
Execution of this load module or phase is terminated. 
System Action With Extended Error Message Fa- 
cility: Execution continues contingent upon the 
error occurrence counts kept in the Option Table. 


7 
cosine, result = 3 (cos x + isIN x). 


IHC2901 GAMMA ARG = X.X (HEX = X), LE 2*** —252 OR 
GE 57.5744 
OA290I 
| ILF290! 


Explanation: In the xxxscAMA subprogram for the 
gamma function, a value of x S 2-75? or x = 57.5744 
is an error. 


Standard corrective action is taken before execution 
continues and consists of setting the functional value 
to *. 
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Appendix D. Storage Estimates 


Appendix D contains decimal storage estimates 
(in bytes) for the library subprograms. The 
estimate given does not include any additional 
library subprograms or FORTRAN execution-time 
routines that the subprogram needs during execu- 
tion. The names of any additional library subpro- 
grams needed are given in Table 15 in the 
column headed “Additional Subprograms.” 

Some library subprograms also require execu- 
tion-time routines for input/output, interruption, 
and error procedures. 

If the programmer has not made allowances 
for the storage required by any of these addi- 
tional routines (see Tables 17, 18 and 18.1), the 
amount of available storage may be exceeded 
and execution cannot begin. The programmer 
must add the estimates for all subprograms and 
routines needed to determine the amount of stor- 
age required. 


System/360 Operating System 


Note: The names of execution-time routines 
sometimes vary according to whether or not the 
Extended Error Message facility is in effect. In 
the following discussion, the names that are used 
are those when the facility has not been specified. 
Table 19 presents a cross listing of names for 
both circumstances. 

The racriosH routine performs input/output 
procedures for both ForTRAN Iv (E) and FORTRAN 
Iv. [This routine refers to a table (macuaTBL) for 
information about the input/output devices used 
during execution.] The rmacrcoMeE routine per- 
forms interruption and error procedures for 
FORTRAN Iv (£) library subprograms; the 
IHCFCOMH, IHCFCVTH, IHCFINTH, IHCTRCH, and 
IHCUOPT routines perform the procedures for 
FORTRAN IV library subprograms. If a system 
contains both compilers, the 19cFCOMH-IHCFCVTH 
routines are used. Tables 15 and 16 indicate 
which library subprograms require these execu- 
tion-time routines, 

In addition, several other execution-time rou- 
tines may be needed to resolve external refer- 
ences in a FORTRAN Iv object module. 


1, If a source module specifies direct access 
input/output operations, the compiler gener- 
ates a call to the mHcpIOsE routine. 
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2. At the point that errors are encouraged during com- 


pilation, the compiler generates a call to an error 


Subprogram 
Name 


xxxCLABS 
xxxCLAS 
xxxCLEXP 
xxxCLLOG 


xxxCLSCN 
xxxCLSQT 
xxxCSABS 
xxxCSAS 
xxxCSEXP 
xxxCSLOG 


xxxCSSCN 
xxxCSSQT 
JHCFAINT 
xxx FCDXI 
xxx FCXPI 
xxxFDXPD 
xxxF DXPI 
IHCFIFIX 
XxxFIXPI 
xxx FMAXD 
xxx FMAXI 
xxx FMAXR 
IHCFMODR 
IHCFMODI 
xxxF RXPI 
xxx FRXPR 
xxx LASCN 
IHCLATAN 
xxx LATN2 
xxx LERF 
xxxLEXP 
xxx LGAMA 
xxx LLOG 
xxx LSCN 
xxx LSCNH 
xxxLSQRT 
xxxLTANH 
xxx LTNCT 
xxxSASCN 
IHCSATAN 
xxxSATN2 
xxxSERF 
xxxSEXP 
xxxSGAMA 
xxxSLOG 
xxxSSCN 
xxxSSCNH 
xxxSSQRT 
xxxSTANH 
xxxSTNCT 


Decimal Estimates 


e@ Table 15. Mathematical Subprogram Storage Estimates 


Additional 
Subprograms 


xxxLSORT 

xxxLEXP, xxxLSCN 
xxxCLABS, xxxLSORT, 
xxxLLOG, xxxLATN2 
xxxLEXP, xxxLSCN 
xxxLSQRT 

xxxSSQRT 

xxxSE XP, xxxSSCN 
xxxSABS, xxxSSQRT 
xxxS LOG, xxxSATN2 
xxxSEXP, xxxSSCN 
xxxSSQRT 

xxxCLAS 


xxxCSAS 
xxxLLOG, xxxLEXP 


xxxSLOG, xxxSEXP 
xxxLSQRT 
xxxLEXP 


xxxLLOG, xxxLEXP 


xxxLEXP 
xxx LEXP 


xxxSSQRT 


xxxSEXP 


xxxSLOG, xxxSEXP 


xxxSEXP 


xxxSEXP 





routine (IHCIBERR for FORTRAN Iv (E) and IHCIBERH 
for FORTRAN Iv). If execution of the load module 
is attempted, the error routine is called, a message 
is issued, and the execution is terminated. 


3. If a FORTRAN Iv (£) source module contains a com- 
puted co To, the compiler generates a call to the 
HCCGOTO routine. 


4. If a FORTRAN Iv source module contains any input/ 
output operations that refer to a NAMELIST name, 
compiler generates a call to the ICNAMEL routine. 


5. If a FORTRAN Iv source module uses the debug 
facility, the compiler generates a call to the nscpBUG 
routine. 


6. If boundary alignment was specified during system 
generation, the rncADJsT routine will be loaded if a 
boundary-alignment error occurs. 


Model 44 Programming System 


In the ForTRAN library of the Model 44 Programming 
System, the BoaFIocs routine is the interface with the 
system input/output services. This routine refers to a 
table, Boaunrrs, for information about the input/output 
devices used during execution. The BoarBcoM routine 
performs interruption and error procedures. If a source 
program contains any input/output operation(s) re- 
ferring to a NAMELIST name, the compiler generates a 
call to the BOANAMEL routine. 


Disk Operating System 


The ILFFI0cs routine performs input/output proced- 
ures for FORTRAN Iv under the Disk Operating System. 
This routine refers to a table, either mFUNTAB or 
ILFGHTAB, for information about the input/output 
devices used during execution. The mFrscom, 
ILFADCOM, and ILFFINT routines perform input/output 
preparation, conversion, and interrupt processing for 
the library subprograms. 

The mLFACoM routine is used, when necessary, for 
communication with pos Basic FORTRAN modules. 
ILFTRBK provides traceback diagnostic messages. 

In addition, several other execution time routines 
may be needed to resolve external references. 


1. If a FoRTRAN module requires direct access input/ 
output operations, the compiler generates a call to 
ILFDIOCS. 

2. When errors are encountered during compilation, 
the compiler generates a call to 1LFIBERR. If execu- 
tion of an incorrect statement is attempted, this 
routine is called, a message is issued, and execution 
is terminated. 

3. If a FORTRAN source module specifies a NAMELIST 
operation, the compiler generates a call ILFNAMEL 
routine. 

4, If a FORTRAN source module uses the debug facility, 
the compiler generates a call to the ILFDEBUG 
routine. 


® Table 16. Service Subprogram Storage Estimates 


Decimal Estimates | Uses Input/Output 


Subprogram 


and Interruption 
44PS Routines 
80 


xxx F DVCH 
xxx F DUMP 
xxxF EXIT 
xxx FOVER 
xxxF SLIT 





Table 17. Execution-Time Routine Storage Estimates, Oper- 
ating System 


Routine Decimal 
Name Estimate Used By 












































THCADJST 1,156 FORTRAN IV 
IHCCGOTO 60 FORTRAN IV (E) 
IHCDBUG 2,152 FORTRAN IV 









































IHCDIOSE 2,688 (See Note 1) Both 
IHCFCOME 6,196 FORTRAN IV (E 
THCFCOMH 4,168 FORTRAN IV 
JHCCOMH2 520 FORTRAN IV 
JHCFCVTH 4,688 FORTRAN IV 
JHCFIOSH 3,744 + IHCUATBL | Both 

(See Notes 2 and 3) 
IHCIBERH 294 FORTRAN IV 
IHCIBERR 136 FORTRAN IV (E) 
IHCNAMEL 2,880 FORTRAN IV 
JHCTRCH 792 FORTRAN IV 


IHCUOPT 8 FORTRAN IV 


Note 1: This module also acquires dynamic storage. Its 
amount, in bytes, may be computed by the formula 
184 + buffer size 
Each buffer value should be 800, except for a card 
read punch which is 80 and a printer which is 133 
FORTRAN utilizes double buffering. 
Note 2: This module also acquires dynamic storage. Its 
amount, in bytes, may be computed by the formula 
128 + buffer size 
Buffer lengths are listed in Note 1. 
Norte 3: The number of bytes in table IHCUATBL may be 
computed by the formula 
12n +8 
where n is the number of data set reference numbers 
requested during system generation. 
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Table 19. Execution-Time Routines Storage Estimates With 
Extended Error Message Facility, Operating System 


@ Table 18. Execution-Time Routine Storage Estimate, Model 
44 System 


BOADIOCS 688 


Corresponding 
Module With No 
Routine Error Message 


BOAFIOCS 1,776 
BOAIBCOM 10,432 
BOANAMEL 2,432 
BOAUNITB (see Note 1) 
BOAUOPT 8 
BNXADJST 1,000 


(Note 2) 


Note 1: The number of bytes in CSECT BOAUNITB may 
be computed by the formula 
8n +8 
where n is the number of data set reference numbers 
requested during system construction. 
BNXADJST does not reside in the library of relo- 
catable modules, but in the absolute phase library. 





Name 


IHCADJST 
IHCDBUG 
IHCECOMH 
THCCOMH2 
IHCEDIOS 
IHCEFIOS 


JHCEFNTH 
ITHCERRM 
IHCETRCH 
IHCFCVTH 
IHCFOPT 
IHCIBERH 
IHCNAMEL 
THCUOPT 


Facility 


IHCADJST 
IHCDBUG 
IHCFCOMH 
IHCCOMH2 
IHCDIOSE 
IHCFIOSH 


IHCFINTH 


IHCTRCH 
IHCFCVTH 


IHCIBERH 
IHCNAMEL 
IHCUOPT 


Storage Estimate 


1,156 
2,152 
5,368 
1,120 
3,848 (See Note 1) 
4,584 + IHCUATBL 
(See Notes 2 and 3) 
1,368 
1,512 
706 
4,688 
824 
224 
2,880 
800 + 8n (See Note 4) 


@ Table 18.1. Execution-Time Storage Estimates, Disk Operating Note 1: This module also acquires dynamic storage. Its 
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System 


Routine Decimal 
Name Estimate 


ILFACOM 
ILFADCON 
ILFDEBUG 
ILFDIOCS 
ILFFINT 
ILFFIOCS 
ILFGHTAB 
ILFIBCOM 
ILFIBERR 
ILFNAMEL 
ILFTRBK 
ILFUNTAB 


amount, in bytes, may be computed by the formula 
184 + buffer size 
Each buffer value should be 800, except for a card 
read punch which is 80 and a printer which is 133. 
FORTRAN utilizes double buffering. 
This module also acquires dynamic storage. Its 
amount, in bytes, may be computed by the formula 
128 + buffer size 
Buffer lengths are listed in Note 1. 
The number of bytes in table IHCUATBL may be 
computed by the formula 
12n + 8 
where n is the number of data set reference num- 
bers requested during system generation. 
The number of additional entries supplied in the 
Option Table during system generation is repre- 
sented by n. 





Appendix E. Assembler Language Information 


The mathematical and service subprograms in the 
FORTRAN Iv library are available to the assembler lan- 
guage programmer. The following text explains the 
method of calling a library subprogram in an assembler 
language program, and then gives additional informa- 
tion necessary to use each type of subprogram. (The 
assembler language programmer should also be fami- 
liar with the information contained in Appendix D.) 


Calling Sequences 


To call either type of library subprogram, the assem- 
bler language programmer supplies an entry name, 
an argument list, and an area used by the subprogram 
to store information (i.e., a save area). The following 
conventions must be observed when calling a library 
subprogram in an assembler language program: 

1. The address of the entry name must be in general 
register 15. 

2. The address of the point of return to the calling 
program must be in general register 14. 

3. The address of the argument list must be in general 
register 1. 

4, The argument list must be assembled on a full-word 
boundary; it consists of one 4-byte address constant 
for each argument. The last argument must have 
a 1 in its high order bit. 

5. The address of the save area must be in general 
register 13. 

6. The save area must be assembled on a full-word 
boundary. Although the minimum size of the save 
area depends upon the subprogram, the program- 
mer is advised to use a save area of 18 full-words 
for all library subprograms. The minimum save 
area sizes are given in Tables 2 through 6 for the 
mathematical subprograms, and in Table 17 for the 
service subprograms. 

7. If the information in a floating-point register is to 
be retained, the programmer must save and restore 
the contents of the register. The subprograms that 
make use of the floating-point registers contain no 
provisions for saving the information. 

8. If a main program in assembler language contains 
any calls to those library subprograms that use the 
FORTRAN execution-time routines (see Appendix D), 
the following instructions must be included before 
the call to the subprogram is issued: 









L 15,=V(IBCOM#) 
BAL 14,64(15) 


EXTRN IBCOM# 
L _-15,=A(IBCOM#) 
BAL  14,64(15) 


These instructions cause the initialization of return 
coding and the interruption exceptions described in 
Appendix C. If these instructions are omitted, the 
occurrence of an interruption or an error causes un- 
predictable termination of the execution of this load 
module. 


Nore: In an assembler language program, a decimal- 
divide exception may occur. This causes the char- 
acter B to appear in the program interruption 
message described in Appendix C. 


The user of System/360 Operating System may use 
several methods to call a FoRTRAN library subprogram: 
the appropriate macro-instructions described in the 
publication IBM System/360 Operating System: Su- 
pervisor and Data Management Macro-Instructions, 
Form C28-6647 or the general assembler language 
calling sequence (given in Figure 3). If the macro- 
instructions are used, the address of the save area must 
be placed in general register 13 before using a macro- 
instruction to give control to the subprogram. For 
example, if the square root of the value in AMNT is to 
be computed and save is the address of the same area, 
the following statements could be included in an as- 
sembler language program to call the micssgrr sub- 
program: 


L 15,=V(IBCOM#) 

BAL 14,64(15) 

LA "  13,.SAVE 

CALL SORT,(AMNT),VL 
inate 


SAVE DS 18F 
If the general assembler language calling sequence 
shown in Figure 3 is used, the programmer must en- 
sure that all of the conventions discussed previously 
are followed. For example, to call the micsserT sub- 
program to compute the square root of the number in 
AMNT, the following statements would be included in 

the source program: . 


LA 13,SAVE 

LA 1,ARG 

L 15,ENTRY 

BALR 14.15 
ENTRY DC V(SQRT) 
SAVE DS  18F 
ARG - DC X‘80’ 

DC AL3(AMNT) 
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arglist 


* 


CNOP 
DC 


* 


13, area 
I, arglist 
15, entry 
14, 15 


Xid’ 
* 


V (entry name) 


A (entry name) 
* * 


xxF 


X’80’ 


General register 13 contains the address of the save area. 
General register 1 contains the address of the argument list. 
General register 15 contains the address of the subprogram. 


General register 14 contains the address of the point of return to the 
calling program. 


This statement is optional. The id represents the binary calling sequence 
identifier. This number is supplied by the programmer and may be any 
hexadecimal integer less than 216 — 1. 


Note: In this case, the entry name must be defined by an EXTRN instruc- 
tion to obtain proper linkage. 


This statement defines the save area needed by the subprogram. The xx 
represents the minimum size of the save area required; however, the pro- 
grammer is advised to use a save area of 18 full-words for all subprograms. 
(The minimum save area requirements are given in Tables 2 through 6 
for the mathematical subprograms and in Table 16 for the service sub- 
programs. ) 


Aligns the argument list at a full-word boundary. 
Indicates the first byte of the only argument. 


DC AL3 (arg) 
or for more than one argument: 
arglist DC A (arg:) 


DC A ( arge ) 


DC X’80” 
DC ALS (argn) 


Figure 3, General Assembler Language Calling Sequence 


When the load module is executed, the micsgrr sub- 
program is called to compute the square root of the 
number in AMNT; the result is stored in floating-point 
register 0. The binary calling sequence identifier is 
not used. 

The assembler language user of the Model 44 Pro- 
gramming System will use the calling sequence given 
in Figure 3. He will, however, use only the A-type 
address constant where the choice between that and 
the V-type is given. As the note in the figure states, 
the label in the operand portion of the address constant 
must be made the object of an ExTRN statement to 
obtain proper linkage. 


Mathematical Subprograms 


The assembler language programmer supplies one or 
more arguments for each mathematical subprogram. 
The arguments may be either integer values or nor- 
malized floating-point real or complex values. 

An integer argument is four bytes in length and 
starts on a full-word boundary. A real argument is 
either four or eight bytes in length. The four-byte 
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Contains the address of the argument. 


Contains the address of the first argument. 


Contains the address of the second argument. 


Indicates the first byte of the last argument. 


Contains the address of the last argument. 





argument starts on a full-word boundary. The eight- 
byte argument starts on a double-word boundary and 
occupies two adjacent words. The first word contains 
the most significant digits. This word is also the ad- 
dress of the entire argument; the second word con- 
tains the least significant digits. 

A complex argument is either eight or sixteen bytes 
in length and starts on a double-word boundary. The 
first half of the argument contains the real part of the 
complex argument; the second half contains the 
imaginary part. The address of the real part of the 
argument is the address of the entire argument. 


Each mathematical subprogram returns a single 
answer. This answer is either an integer value or a 
normalized floating-point real or complex value. An 
integer answer is stored in general register 0, a real 
answer is stored in floating-point register 0, and a 


complex answer is stored in floating-point registers 
0 and 2. 


Tables 2 through 6 contain additional information for 
using the mathematical subprograms in an assembler 
language program. These tables give the floating-point 


Table 20. Assembler Information for the Service Subprograms 


Subprogram Entry Save Area 
Name Name(s) (Full — 


xxxF DUMP DUMP 



















PDUMP 
xxx FDVCH DVCHK 
xxxF EXIT EXIT 
xxxFOVER OVERFL 
xxxF SLIT SLITE 





SLITET 





registers that are used by the subprogram and the 
save area required by the subprogram. 


Service Subprograms 


The service subprograms do not use the floating-point 
registers during execution; however, each service sub- 
program requires a save area. The minimum size of 
the save area depends upon the subprogram to be 
used and is given in Table 20. 
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Appendix F. Sample Storage Printouts 


A sample printout is given below for each dump for- 
mat that can be specified for the xxxrpumMP subpro- 
gram. The printouts are given in the following order: 
hexadecimal, LOGICAL *1, LOGICAL *4, INTEGER *2, 
INTEGER *4, REAL *4, REAL *8, COMPLEX *8, COMPLEX *16, 
and literal (see Figure 4). Note that the headings on 
the printouts are not generated by the system, but 
were obtained by using FORMAT statements. 


The output of the xxxFpumMpP subprogram (for both 
entry names, DUMP and ppuMP) is placed on the object 
error unit data set defined by the installation during 
system generation. 


CALL PDUMP WITH HEXADECIMAL FORMAT SPECIFIED 
OOAZEO 485FS5E10 00000000 485F5E10 10000000 42200000 


006D0C8 42800000 00000000 00000000 00000000 00000000 00000000 00000000 00000000 o0000000 0000000 0000000 o0000000 
0O06DF8 C0000000 90000000 41200000 41566666 od000000C 41100000 


CALL PDUMP WITH LOGICAL*1 FORMAT SPECIFIED 
OO6E1E T OF 
CALL PDUMP WITH LOGEICAL*4& FORMAT SPECIFIED 
006E10 F oT 


CALL POUMP WITH INTEGER*2 FORMAT SPECIFIED 
0O06E18 10 

OOGEIA 

OO6EIC 

CALL PDUMP WITH INTEGER*4 FORMAT SPECIFIED 


006E20 
006E48 


CALL PDUMP WITH REAL*4& FORMAT SPECIFIED 
Oo6E0O 0.20000000E 01 0.53999996E 01 
CALL PDUMP WITH REAL*8 FORMAT SPECIFIED 
006DC8 . 0.1759999999999999D 03 

CALL PDUMP WITH COMPLEX*8 FORMAT SPECIFIED 


006DD0 €3.0000000,4,.0000000) C€4,0000000,8,0000000) 


" CALL PDUMP WITH COMPLEX*16 FORMAT SPECIFIED 


006DE0 €0.9999999999999990, 0.9999999999999990) C€-0.9999999999999990,-0,9999999999999990) 
CALL PDUMP WITH LITERAL FORMAT SPECIFIED 


QO6ESC THIS ARRAY CONTAINS ALPHAMERIC DATA 





Figure 4. Sample Storage Printouts 
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Absolute error ........... 00000 e eeu Se eens 21, 41 
Absolute value) 1. oth eawiddasnee teres etee es 7, 11, 22, 41 
Accuracy Statistics: 2250562 p04 4.18 eee wea eme wed 41, 47 


AINT (see IHCFAINT ) 

ALGAMA (see xxxSGAMA) 

AVGOPEBIMNS 4.5: oho eh eae se 2 ha Peewee Nees 21, 40 
ALOG (see xxxSLOG) 

ALOG1O (see xxxSLOG) 

AMAXO (see xxxFMAXI) 

AMAXII (see xxxFMAXR) 

AMINO (xxxFMAXI) 

AMINI (see xxxFMAXR) 

AMOD (see IHCFMODR) 


Arccosine subprograms .................-. 7, 9, 24, 32-33, 56 
ARCOS (see xxxSASCN ) 

ATOUMENES 3 io we odes dhs RGR CAE eT Rs Ga Re ROE 6, 60 
ARSIN (see xxxSASCN ) 

Arcsin subprograms....................5- 7, 9, 24, 32-33, 56 
Arctangent subprograms............... 7, 9, 25-26, 33-34, 56 
Assembler language calling sequence................. 59-61 
Assembler requirements .................. 7, 8-13, 14, 59-61 


ATAN (see IHCSATAN or xxxSATN2) 
ATAN2 (see xxxSATN2) 
CABS (see xxxCSABS ) 


Calling sequence ..........0 0000. ccc eee 59-60 
Calling FORTRAN subprograms 
SXPUCHY 6 o Pogs a aye eGrle awe hein Renda wee 6-13 
THIPICILY 6 2 ep eacrce 2 Qe dns Phy aloe x aairlneg eae 14-15 
in assembler language .................0 0002005 59-60 
CALL macro-instruction ................2.-..-00005 59 
GAIL Statement oie 33 oot cores ade ow et ed-ebd 5, 16-17 


CCOS (see xxxCSSCN ) 

CDABS (see xxxCLABS ) 

CDCOS (see xxxCLSCN ) 

CDDVD# (see xxxCLAS) 

CDEXP (see xxxCLEXP) 

CDLOG (see xxxCLLOG) 

CDMPY# (see xxxCLAS ) 

CDVD (see xxxCSAS ) 

CDSIN (see xxxCLSCN ) 

CDSQRT (see xxxCLSQT) 

CEXP (see xxxCSEXP ) 

CLOG (see xxxCSLOG ) 

CMPY# (see xxxCSAS ) 

CSQRT (see xxxCSSQT) 

Common logarithm subprograms... .. . 7, 8, 21, 28-29, 36-37, 56 

Complemented error function subprogram..... 7, 11-12, 26-27, 
34-35, 56 

Corrective action after program interrupt occurrence. . . 18 

COS (see xxxSSCN) 

COSH (see xxxSSCNH ) 


Cosine subprograms.......... 7, 9-10, 23-24, 29-30, 37-38, 56 
COTAN (see xxxSTNCT) 
Cotangent subprograms............... 7, 10, 31-32, 38-39, 56 


CSIN (see xxxCSSCN ) 
DARSIN (see xxxLASCN ) 
DARCOS (see xxxLASCN) 
DATAN (see IHCLATAN or xxxLATN2) 
DATAN2 (see xxxLATN2) 
DCOS (see xxxLSCN ) 
DSCOSH (see xxxLSCNH ) 
DCOTAN (see xxxLTNCT) 
DERF (see xxxLERF ) 
DERFC (see xxxLERF ) 
DEXP (see xxxLEXP) 


Index 


DGAMMA (see xxxLGAMA) 
Disk Operating System 


Assembler Language .............- 0.000 e eee ees 2 
Error Messages ....... 00.00 c cece ees 50-55 
Execution-Time Routines ................0000005 57 
Interrupt Procedures .........0.. 000000 eee eee ee 49 
Storage Estimates ..........0. 000.00 c cee eee 56-58 
Subprogram Names .............. 0.00000 e ee eeee 5 
Divide-check exception .......... 0.000000 cee eens 16, 48 


DLGAMA (see xxxLGAMA) 
DLOG (see xxxLLOG) 
DLOGI1O (see xxxLLOG) 
DMAXI (see xxxFMAXD ) 
DMINI (see xxxFMAXD ) 
DMOD (see IHCFMODR) 
DSIN (see xxxLSCN ) 
DSINH (see xxxLSCNH ) 
DSQRT (see xxxLSQRT ) 
DTAN (see xxxLTNCT) 
DTANH (see xxxLTANH) 
DUMP (see xxxFDUMP) 
DVCHK (see xxxFDVCH ) 
Entry AME oor ee a EASES oe OS BRS 6 
ERF (see xxxSERF ) 

ERFC (see xxxSERF ) 


Error 
absolute .......... 00. ee eens 91, 41 
MMESSAMES 5 ik Me wees Se ae ease phate et 49-55 
Optional Service: hb s ga 44 eee be us tg Shit ee ei 16 
procedures. .cah) hy se eG Phos hee ee ee as 49-55 
DFOPAGAUON, fyi. bose ds Oe te we ek Re eee ees Al 
POLATIV GY 6 ooh eae ee Aa eles Pa 21, 41 

Error function subprograms......... 7, 11-12, 26-27, 34-35, 56 

Execution-time routines ...............00-00 020s eee 56-58 


EXIT (see xxxFEXIT ) 
EXP (see xxxSEXP ) 


Explicitly called subprograms ...................04- 5, 6-13 
MIS OF: epee oAe-c Phe bs tie Be tank AE ROR CREAR EES 7 
performance statistics .......... 0.000. c eee ee eee 42-46 
SIZE: “OF bys fhm Alert eae DESY ese ae ERT Ba 56-57 
fabless coc AN Se dace iid te MAG ta a A ie ee 8-13 
sein: FORTRAN. (35. csvesreeae ee ae aR E SS 6-7 
use in assembler language.............-...-..-.-- 59-61 

Exponential subprograms ................ 7, 8, 22, 27, 35, 56 

Exponent overflow exception ................000005- 19, 48 

Exponent underflow exception ...................004. 19, 48 

Extended Error Message facility..................... 15 


FCDXI# (see xxxFCDXI) 
FCXPI# (see xxxFCXPI) 
FDXPD# (see xxxFDXPD) 
FDXPI# (see xxxFDXPI) 
FIXPI# (see xxxFIXPI) 
FRXPI# (see xxxFRXPI) 
FRXPR# (see xxxFPXPR) 


FameGon “Vale: 6.5 sedis eves PORE ON Oe na ee as 5,6 

GAMMA (see xxxSGAMA ) 

Gamma subprograms ...............++0005- 7, 12, 28, 36, 56 

Hyperbolic cosine subprograms ............. 7, 11, 30, 38, 56 

Hyperbolic sine subprograms ............... 7, 11, 30, 38, 56 

Hyperbolic tangent subprograms ............ 7, 11, 31, 39, 56 

IDINT (see IHCFIFIX ) 

Implicitly called subprograms ................... 5, 6, 14-15 
TSE ihe cutie ik hc ta ae Gera RNS ate EATS ocho Og ete, dese 14 
PESule OF Use ieee aa ee OE ER 15 
SIZE ita eel ee Ae eG NE ig ta Ratan! et anche adie 4, 56 
USE) i eet ee oe etn ie ates sd Sashes ee ees 14 


INT (see IHCFIFIX) 


Interruption procedures ............. 0.0000 c eevee AT 
Linkage-€ditor oe cos. eae eS ee eee ee REE tees a8 5 
Logarithmic subprograms........... 7, 8, 22, 28-29, 36-37, 56 
Log-gamma subprograms .................0000005- 28, 36, 56 
Machine indicator test subprograms............... 19, 58, 60 
Mathematical subprograms ....................-4-. 5,6 
BICOMLAINS ho oho aeate Roa a Rat eee 21, 40 
CEHMIIONS. 2% 2h tho-5 aes haces eee asaa Ep Sale’ 5 
explicitly Galled \.).0.3. Saaeiee aba es ated ketenes 6-13 
implicitly Called: c:c2. facsu suis dota een hee es 14-15 
MSE IOP. 22s 6 Ao Slut a neat baba oe Bo ey eee anew he 8 7,14 
performance wesc enh Ses os fae oe ee eee e RS 41-46 
SIZES Sion eat ence Ae ockdes Ue arate he eee ws 56 
use in FORTRAN 0 4.66:5 bce cen kee ah eeea eden 6-15 
use in assembler language....................... 60 
Maximum value subprograms ................... 7, 12-13, 56 


MAXO (see xxxFMAXI) 
MAX] (see xxxFMAXR) 
MINO (see xxxFMAXTI) 

MINI (see xxxFMAXR) 


Minimum value subprograms ................... 7, 12-13, 56 
MOD (see IHCFMODI) 

Model 44 Programming System..... 5, 41-46, 47, 49, 58, 59, 60 
Modular arithmetic subprograms ................... 7, 13, 56 
Natural logarithm subprograms .................... 7, 8, 22 
Operating System, System/360......... 5, 41-46, 48-49, 56, 59 


OVERFL (see xxxFOVER) 
PDUMP (see xxxFDUMP) 


Program interrupt corrective action.................. 18 
Pseudo sense lights ........... 00.00. c cece eee ences 19 
Relative: Gor 2s ssck ate devesciacean teu cf ese veees 21, 41 
Sample dump printouts .................0 000.0000. 62 
Sampling techniques ............0....0 0000 eee ee eee Al 
DEUSEISHES Ac caer ge ature lttee elated Anitesh eecald a 19 
Service subprograms 
machine indicator test ................0 00200000. 19 
SIZES) x svcscelereits Lane ety ae ee hat he WN ORE A 57 
use in assembler language....................... 59-62 
psec MORDTRAN a 6 dc tt sa, tiple edn wate es os 19-20 
Wilby teehee weed ee nae eho des cea 
SIN (see xxxSSCN ) 
Sine subprograms ............ 7, 9-10, 23-24, 29-30, 37-38, 56 


SINH (see xxxSSCNH ) 
SLITE (see xxxFSLIT) 
SLITET (see xxxFSLIT ) 
SORT (see xxxSSQRT ) 


Square root subprograms ........... 7, 8, 22-23, 30, 37-38, 56 
Standard deviation .........00..00 00 0c ccc cece cease Al 
Storage estimates ........0 000.00. c cece ec eee eee 56-58 
Storage printouts ....... 20.0.0... ccc cece cece cee neuas 62 
Subprogram names ............. 000000 c cece eeueee 5 
Subprograms and execution-time routines 
BOADIOCS routine ................0 0.000000 e ee 58 
BOAFIOCS routine ................ 0.0.00 eee. 58 
BOAIBCOM routine ................ 000 0ce eee, 58 
BOANAMEL routine ......................000.. 58 
BOAUNITB routine ............ 0.0.0... ..00005. 58 
BOAUOPT routine ............. 0.0.0.0 0 cece eee 58 
BNXADIST routine ................... 0000-000. 58 
IHCCGOTO routine ................... 0000005. 57 
xxxCLABS subprogram 
BIGOMUNMY hg. 2, th cudnt, Soe cant dh heh aemumee ewes: 22 
effect of an argument error .................... 22, 
performance ............. 000 cece cece ee eee 42 
SIG? Sc Bio th eat he Seana Sixes wera aeaieces Ua edeae 56 
USO 8 2: te tosh DiS Sas ce acta ues ene tee ui eet A Pm ees Bt 11 
xxxCLAS subprogram 
SIZE ids see eect ee fom Dak Pilea died boats hee 56 
USC Fi tite attic wtb eea wl sray ane eines 14 


64 


xxxCLEXP subprogram 
AloOrithaa: iS sa Real eRe Eee we Blan, oats tha tis 27 
effect of an argument error .............. 0.00005 27 
GITOT ANGSSAVES. | oN whe Gnd Raabe haw Bea 52, 54 
performance .......0 0.0.0 cece eee eee 42, 
S120 2 oa ia UB Aeas Peet pA ae unten coeds 56 
MSO mere xsi a subd Ba tadte eee Gt A ches oe tat B oe, Lick ates cane 8 
xxxCLLOG subprogram 
I POTIEINL 25.225 avs £4 Aoi caecawa ge Sp ne erage ene eins 22 
effect of an argument error..................4.4. 22, 
CITOR MESSAGE 5.05 va yok Ra he bs So week RRS ORS eR 54 
PETIOLDIGNCE epee os cae hea ee eS A2 
SIAC! adns aig oa eta BSA od ke Mp eRe Oe ee ee 56 
VSO. athe te ah a AP ocd ans Dien ies Ge he hee ad 8 
xxxCLSQT subprogram 
alPOMtAM, 54M aaa atoedurg ok he mua ae dx 22-23 
effect of an argument error..................00- 23 
PerormMance. ai yah. cule wad GR Se A we eae a ey 42 
SI ZO isos ks eit, Ody ake eran? eeiee andy de AC aloe 56 
WSO? snd on TAD HF Doin hel i EER ee ot ER es 8 
xxxCLSCN subprogram 
BROLIN an ep OR RG eh Sr ee ee wea 23 
effect of an argument error..................... 23 
CLVOD MESSAGES eh ol Sos ET he Mae hee See ty ek 51, 54 
PSYFOMMANCE 1 secs acteysayin's te Glo eee eae Saeed wee ai oe 42, 
SIZE! 4664. Saad Sota Bees tebe a I, aati 56 
MISC eee Sr et uch ads ho ce Bhd o Aas Inlnte aioe Gun EA os 9 
xxxCSABS subprogram 
MOTI: oSicaty. ce etch acs de Nan some hat pagel tens. a asks 22 
effect of an argument error..................04. 22, 
PerlOMnAnee- le le 40s cat ee ake ee eee 42 
SU ps hel acedeieai Miah aa a ethie Pree Nek S aateaee cath wid 56 
MISC hes snus Sek oni a aye ete SMa tak Rea ate wen ll 
xxxCSAS subprogram 
SZC? hiked Gane tht eee ee teas le bt cay Gea tae Bo 56 
USC sate Shaler aya xi Rang cesta sent ad ino se Tota ala ean 14 
xxxCSEXP subprogram 
BIG OMION: ini wine BO re eaten oh Gkndled Atel Be hceende 22 
effect of an argument error..................... 22 
CITOL MESSAZES ys. e en edn ae ee Seo eee ESA ESS 51, 53 
DETLOTIMANCE: 546 bid 5 8G Sak dd Leah eee ele has 42 
SIS: ag Si tee er ese dee Seer bne gs Raed Salas Sond 56 
USE So fa tienes ead fd BRR OS Ue Se ed 8 
xxxCSLOG subprogram 
BIB ORIRNT ad oo tare itless @ Oia de tues y Ym eke ea ky 22 
effect of an argument error..................... 22 
GITOR-INESSACES Sb Siss. 2.5 oe ee ONG bu ee eae 53 
DEVLOLINANCE. £.o un he-eS ke oe ees ee ke Rob Sak 42 
SIZ is ecta tiie A tos occ, oh a A cee tse Aya cade eee 56 
MSOS 658g BREE Le ht O54 bo oe oe dee Sea’ 8 
xxxCSSQT subprogram 
SIGORILTEe iea i. oe ig, saidis idee alan ele as eeee 223-23 
effect of an argument error..................0.. 23 
performance: 45:4 h5. 6-0 05 oN eee a dds 43 
SIZE if acco Gx dw da Gy ahah NE aed ead fa ee phe 56 
NIGGA sie leony arate < rout ento Monaitarliad aa oe ead es aeacs ese ai 8 
xxxCSSCN subprogram 
BIOTIN 6a ole op aes ena beeen eek heme bette? 23-24 
effect of an argument error..................... 24 
CIMPOr MESSAGES 66 s83i5 sa ee eek ee eeu dg abicur aaah 53 
performance ............ 0.0... cece eee ee eae 42-43 
SIZE e255 ad, eit ls ag ate ne Paced ete eee 56 
MSO? iain Paha oe Ste Mara hander hate: sehaiandcondl eas ate aul: othieths 9 
IHCDBUG routine’. 9510.40 Sars o eae cee ants 56-58 
IHCDIOSE routine ............. 00.0.0 .0 cee 56-58 
IHCFAINT subprogram 
BIZ G gsc su Wks aed Sh Sau RG Weal dy ie dda ee SA eae 56 
USC: sone; sens Bisa hae Gee eke CR aie VES ote 13 
xxxFCDXI subprogram 
GITOL “MESSAGC® oe sab hates wo he eR aaa 50 
TOSUIE OL USC: eines ens 7s Kee tenet sure ile awa del 15 


SEZOl ethyl BW een ake hdd ate ek 

BISEE epcicgcahs Dkk te erin eB Ghee int tare De A Amey ase 
JHCFCOME routine ................ 0.0 cee eee 
IHCFCOMH routine ................... 0000s 
JHCFCVTH routine ............... 0.00002 eee 
xxxFCXPI subprogram 

EYTOr MESSAGE .... 6... eee eens 

PESUIE GE USE 6:5. 30c et en a BAD eee RS 

SIZE Paced tens hoods Hie BSG a nk adn AI tse by ake 

BSc in olan Sead ie aE arse caso Baul onde a ENS 
xxxF DUMP subprogram 

assembler requirements ................-..005. 

format specification ............0...-0000 eee eee 

OUtPUE 4 oct ana aati eb Mategg Boo aa dee 

programming considerations ................... 

sample printouts ...............0 0.0 cee eee ees 

SIZ Gs sari i atheche hen Gk susie ASE Geet a MRS Net a fees 

GS: hy cue dexc ahd natin inne eee ahh e Oe ie EE eS 
xxxFDVCH subprogram 

assembler requirements .................0-00005 

SIZE:hs seal sath BSriceel See Babe ORGS 6 age echt 

MISE Bate Saya paMnticety Shaky HIBS bel Salata ceva ye RE 
xxxFDXPD subprogram 

GiTOR MESSARE- bacucka dewalt bd 2 ae aaa dine SA nal 

FOSUMIESOF NSE 1n a ia ee KG nae ak edad 

SIZO 8 ip. h Ue ASA DORAL or PO eh eg a bow ie 

VSG esi RAPE ts ht eh hk eas os Go AOR OG 
xxxFDXPI subprogram 

EITOr MESSAGE 2 oo eee eee ese eeees 

Pesull OF AIS 56.55% shai keg Haetig es Bioule Sak Sa ae 

SIZE Hoses ate etic a eis eee heed ee ns Dares 

USE. aa SA EEO odes eating aged ante 
xxxFEXIT subprogram 

assembler requirements ...................004. 

SIZEO sed Shes Bie pty Re aah PEO sok Sten ners Pn aeaNeye hE kas 

USE: oh hse @ ot Lilie ee anak mee eh de teh ye Pena are 
IHCFIFIX subprogram 

SIZE! sae ad os As bag gat Sed Paes oe ee ea a ek 

RASS es GSO tae ile ties BOL os tos gw Ae og eae 
xxxFIXPI subprogram 

error message ........ 0.0.0.0 ccc cece ee eee eae 

PES OL MSE). 92. 2u5 eo ree honored nw tals 

SIZE oases At Canes NN Oa cn Slash oh ere Goa aeamethe ne nog 

USE osc te SaaS eRe eed oeig sd Eos eh wenn 
xxxF MAXI subprogram 

SIZE noth lA hee osdealens Ue aad LR edd oe haps 

HSE Jolene ie Sah PAV Nanas Ue tetas whe ea ay 
xxxFMAXD subprogram 

SIZ Ais are Ring a ictaneae Daal eM raate Mower Shc ants aoe Sahee donee 

USE ge eee nis ghee Deady oor eee beon eee 
xxxFMAXR subprogram 

SIZE pind eae eae whe Peel a ed twee eo 

MSG 3 eich Svat, ie Pat Wa Saag pt tet ata? wick hoes aden Oe 
IHCFMODI subprogram 

SIZE savin apie settles tee dead Paka oe Dad eee ete Seatac 

NSE osc fee SAMAR ON we caluanen eae aes 
IHCFMODR subprogram 

1 a So OO 

MISE) secon tytsader hye oe cig Aa arated he Mia cates TN ats 
xxxFOVER subprogram 

assembler requirements ....................... 

BIZ cssiysie YG Bee uyh 8 die beach Gene ew aiens Loko Dose My 

ISG ns aestoniSon Vise p ae yore) Baal nae dass Gis eek Mi lo yee 
xxxF RXPR subprogram 


error message 
result of use 
size 
use 


xxxFRXPI subprogram 
CLTOR MESS CE. oo. 6 a6) 5 Se da, Rider c rede 2 acl acsta Bs 
result of use ............... 


56 
12 


56 
13 


56 
13 


SIZE ties tet coh irerat ene we aeeto ee a Ree Sat eG BA Bk ed 56 
MSGS 25 isis Phe Seas ea ed eRe ee OS Ree 14 
xxxF SLIT subprogram 
assembler requirements ...................0.0. 53-54 
CITOT MESSAGE. 0.5.5.4 ieee sale bee eke pan aa es 50 
BIO cu byrne Se SE arate and Bib alstdanie wie ixiedne, Avene Lye 
LT = a a ve 16 
IHCIBERH routine ...................00.00008. 56-58 
IHGIBERR: routine: <is5 22 jase keds udude ane ad 56-58 
xxxLASCON subprogram 
PIPOMUI inet ls os he al mitt hentai daalacke oh 24 
effect of an argument error...................5. 24 
error message ........ 0... eee 52 
PETLOTIMANCE op tise ah ance a tin eee race hrc emans Band 43 
BIR. Nos es utes tsa hag Wihiaek (eee ama boa Ae anon LOR 56 
WSO. Sx Sct phi tich icin habeas pad Marto & aay Rod 9 
IHCLATAN subprogram 
SISOMCHI on. 32 2c SAT oho RON saa ae bed bas 25 
effect of an argument error.................... 25 
PETPOMANCe oie ne yrenese pana Farag asad once d 43 
SZC tele heh auee Re I ee ORR Oe ae 56 
WSC) et sitew A ote t Shik OS vnc cates ate a ce Oe ae oe 9 
xxxLATN2 subprogram 
HIGOMMN aise Katey Diy ate a ee ae 25-26 
effect of an argument error..................... 26 
efor Message ties el esd hte etek 52 
PEOIMaNce: «css neuen, WdevdaeeaG ose 43 
SELEY soccer ee dt i POW at ES De OF Pilsen Bes 56 
MISES ene ther apes oe Dede ed ee ROE Oe Badd are 9 
xxxLERF subprogram 
ALCORN «Shs, a5 F aeccoraeh tree Rete Ratt Ses Vacate 26-27 
effect of an argument error...................0. 27 
pertormancé: 9.4.24 sh hea ae Sede eakas 44 
SIZ eid tie tres Spd tat tes Apa BORE RR: ces ih 56 
MSG). eo Aa Paitin ed pas Wwe Renew ae Be ws ed oa 11 
xxxLEXP subprogram 
al eomthy oe aoe rag dtaewisly hee au ele ae bo clee 27 
effect of an argument error.................... 27 
CPrOr Messave bel est be hoy hoe eR BR 52 
performance ......... 0.000000 c cee cee eee 44 
SUZ oy Goi tineh etexed ea BE doen use Oe a es 56 
WSC) 25.4.uii ah ei oie ee chee a vee ees 8 
xxxLGAMA subprogram 
AIPOM NIN Sign chde 3 = elle See oh a dale oe mecca wes 28 
effect of an argument error..................... 28 
error messages .... 0.0.0... 0c cee eee 55 
performance ....... 0.000. c ccc ccc e eens 44 
SIL ey Ok Scrat ai Rca read WAG Atak ok, kg 56 
WISE ios Sohn ile Pes tudo lek Sanaa avn ng Ors Monegan 5 Sai y Woneie x 12 
xxxLLOG subprogram 
ALGOPIENGY 5.5 ich Sloe wee Vas 6 See POS ee 28-29 
effect of an argument error..................... 29 
Orror message f3 v2 Cty so ek ht ae ch eed os 52 
performance ........... 00. ccc cect eens 44 
ZO int dnd al pean. BA a ee Rea ie ate 56 
VSG! cao fol Sass PER Fae a eae aaa ee 8 
xxxLSCN subprogram 
BIGGIN: 0.2.9.6 7 Roatan od Ay wahe ne hos eee es 29 
effect of an argument error..................-. 30 
GIO? MESSAGE aa whee ads ORs EER RS pa 52 
Performance 224-029 4ss wage yp naedws bean habe oees 42, 43 
SIZE Sn nos Aetna sh a pede Careers, CAR Chetan D 56 
MISE? sed ins nie Sod etalaiak Ray ott aa Bown gd Pete Mena wees Eee 10 
xxxLSCNH subprogram 
AGORA Ae toon sole CaS Loe whos arate me eee 30 
effect of an argument error.................... ‘ 30 
CITOL MESSAGE: vs dg oe Mates iad ckas ESE ev eo 52 
pertormance «675.442. 4b Ss Wed owes eet eeeS 43, 44 
BILE ocr ge i aed he oA RS be bee A ES 56 
MSE™ sh cM ectsnatteaes Moth. ths ee eee Ae er en Shee tes 11 


xxxLSQRT subprogram 
al eOrithin: ia: $34 Sawer ERE Medes SESE 
effect of an argument error..................04- 
EITOr ‘MESSAGE «0... eee eee eens 
PErfOTMANCE <- Lay ee vane siete ce eeeer ede ewan Ses 
SIZO- os Seok ye RE Se al ae ers ee le Re Es 
USE ee ici aha bg ae Rea S eee es 


xxxLTANH subprogram 
SIGOVITNM 435 y-scet wa oe awe TR weal a ediacely ae in 
effect of an argument error...............-.-06- 
Performance «. os. 6 6 fs 64 RAL eK ee ESS 
SIZE: he nats Doe ep aan ete ee NOE Gods ogee ea Ae este 
USO so a8. 55N5, Bes. Saralsec bese etd he eee: Ma eco aris BF ere el ee toy 07 oe Aer) a 


xxxLTNCT subprogram 
SPOT: cGy xin vos Ree hens Gale os Dewees esas 
effect of an argument error..................05- 
Error MeESSageS ...... 0... cece eee ee eee 
performance «..46 65 esd whe seca eee eh evens 
SIZ oe as hives Sei avec fotin let S gele:  ht wea Mate Beta: ova geo aeaaS 
MSG 3 Ord ecg Moh see Case ae Da ee RAG eal Bape e 


JHCNAMEL routine ......... 0.0... ccc eee 


xxxSASCN subprogram 
DLGOTIC NN hy hk be as Ram elem care Niet 2 we dee ncetet 
effect of an argument error.................000- 
CITON Message. hese soo eben Ree ee evo ne ewes 
PSTLOMMance’ 6046s boy ede Ree edare de pee A mons 
SIZC Hii o hs eS ata iene AE bpd eeu es Be tart ea 
USO 2 le storie eh Wau Mite wanae bed a p ia iGtow hee 


IHCSATAN subprogram 
AOPILUNY nic eRe a ae ek Ee ea ee Bas 
effect of an argument error..................00. 
PErOrmMance? i, 5.o<4 kas oe tou k NA ete pets 
size 


i 


xxxSATN2 subprogram 
DIGOMENNG. Aca GucaSd ee bet Ree. cuees gee ewecgaees 
effect of an argument error..................05. 
GMOr MGSSAGE. oh 2s ek ce hd eae es RES 
performance ............ 0. cece eee cece eens 
SIZE 686 soles etd po bid eee Soh pcb ates assayed 


xxxSERF subprogram 
SIGOLIENIN. «0406 atin A Ge oS Ral Seek od 
effect of an argument error.................... 
performance ......... 0.0... cece cee eee ee ees 
size 


xxxSEXP subprogram 
IDOLE be fares 'e sd Mike ayant § ou aly at oas 
effect of an argument error...................-. 
error message ........ 0... c cee eee eee 
performance ......... 0.000. c cece cece ees 
size 


a 
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xxxSGAMA subprogram 


algorithm. v6 46 os Mee a oe CEREALS DAR EE 36 
effect of an argument error............-..-0005- 36 
EITOF MESSABES .. 1... ee ce tee eee 54 
PerfOLMANCe 466d. si segs Tae Nae eK ae de eee 42,45 
SIZE Berea eta Co aaca: mike ears A edieg teas Sen eae MS Bid aE 56 
MISE sth etic ecete a eseed wkele fe es rae A aie tee ae 12 
xxxSLOG subprogram 
AMON. 45.8 KADER SE ee Ie ha Pasa 36-37 
effect of an argument error..................055 37 
@rror MESSAGE oh is os had ease as ee ERS 51 
DEMONMANCe: 6 od50 2 Seana eee SRS es 42 
SIZE? eds hr bu awk De ais wee tne Pha eee re 56 
MSE. alesis tb ied at a Ie oie alte a grating Shyla teraiee Gea ete 8 
xxxSSCN subprogram 
AlBOrIEAA: «ses, o ae kA esl pe as WEN oR ae 37-38 
effect of an argument error..................... 38 
€fror MCSSAVE liad ei ou dia eke ew eet eee Eh eee 51 
WELPOTMARCE «x oiadided one eh ei awa tae See ees 43, 45 
SILO nds gee abe sn ered Seep eS ded lw rete acres ws ae nae Ne ee ed 56 
MISOE cls ee 6 Ware rath ea a ayes ly eet, tod asd ao 10 
xxxSSCNH subprogram 
SeOrthm 3, 40% spare ee oa kaye ee gales 38 
effect of an argument error..................04. 38 
EIrOr MESSAPe ... 1... ee eee 51 
PEMOMMANCE: sa8is Soe UN hh hg Gs eS 43, 45 
SIZ6 2 Pid eed ls LOE ee ARS in wae oe Es 56 
USC 3 eee etece Sneha we a alee a wrens eee Il 
xxxSSQRT subprogram 
P SOTILOIN: ; «A Sirgee of BSG ean ante hog evans Hea eae 38-39 
effect of an argument error.................... 39 
error Message ......... 0.0. cece cee eee tenes 51 
PerlOrimance. — sas Si cistern d Shadi waloh emer ee hae 45 
SING 6G cating & heed ate ae tees Boe ok 4 Sa 56 
MIS) sh boxce bie tlerene wine oasis 0s Looe eae 2 ta 8 
xxxSTANH subprogram 
WiPOritha, oe. 675 eB eS Caren th ae as one 39 
effect of an argument error...................4. 39 
Perlormance:? 64.20.26 554g ots ated eka otnwe ens 46 
SIZ 4 hap ned Beas Cons aoe a bee eee Fe ee 56 
USE og igat yee sancti ee aadaccmanid ek thew dake 1l 
xxxSTNCT subprogram 
Bl eOrithin 6.2m Bes arses Dh eet eS le he nn 39-40 
effect of an argument error..................04. 46 
CITOY WWESSAVES Coc isc a bcd eee Ha Oh ee GS 51, 52 
PPEPCOMMANee: Ho so ae ec oe a oe eee Salter s 43, 46 
SIZE) hae oh he a pe ee ae Rew a eed 56 
MISO Gok SS sh Al a ea OSe eae eas FS 10 
TAN (see xxxSTNCT ) 
Tangent subprograms................. 7, 10, 31-32, 39-40, 56 
TANH (see xxxSTANH ) 
‘Paming statistes: 22 2s Sata vegeta crag eee ands 41-46 


Trigonometric subprograms. . . .7, 9-10, 23-26, 29-34, 37-40, 56 


Truncation subprograms ..................000000- 7, 13, 56 
User-supplied corrective action .................... 16-18 
Utility subprograms ...................055 16-17, 56, 58-61 





C28-6596-4 


TBI 


International Business Machines Corporation 
Data Processing Division 

112 East Post Road, White Plains, N.Y. 10601 
[USA Only] 


IBM World Trade Corporation 


821 United Nations Plaza, New York, New York 10017 
{International} 


READER'S COMMENT FORM 


Title: IBM System/360 Form: C28-6596-4 
FORTRAN IV 
Library Subprograms 


Your comments assist us in improving the usefulness of our publications; they 
are an important part of the input used for technical newsletters and revisions. 


Please do not use this form for technical questions about the system; it only 
delays the response. Instead, direct your technical questions to your local 
IBM representative. 


Corrections or clarifications needed: 


Page Comment 


Please indicate in the space below if you wish a reply: 


Thank you for your cooperation. No postage necessary if mailed in the U.S.A. 


8-6596-4 


‘OUR COMMENTS, PLEASE . 


[his publication is one of a series that serves as a reference source for systems analysts 
yrogrammers, and operators of IBM systems. Your answers to the questions on the back of 
his form, together with your comments, will help us produce better publications for your 
ise. Each reply will be carefully reviewed by the persons responsible for writing and 
oublishing this material. All comments. and suggestions become the property of IBM. 


*>LEASE NOTE: Requests for copies of publications and for assistance in utilizing your IBM sys- 
em should be directed to your IBM representative or to the IBM sales office serving your 
ocality. 


our] ST) Buore yno 


FIRST CLASS 
PERMIT NO. 33504 
NEW YORK, N.Y. 


BUSINESS REPLY MAIL 


NO POSTAGE NECESSARY IF MAILED IN THE UNITEQ STATES 





POSTAGE WILL BE PAID BY... 


IBM CORPORATION 


1271 AVENUE OF THE AMERICAS 
NEW YORK, N.Y. 10020 


\ttention: PUBLICATIONS 


nternational Business Machines Corporation 
Jata Processing Division 

12 East Post Road, White Plains, N. Y. 10601 
USA Only] 


BM World Trade Corporation 
121 United Nations Plaza, New York, New York 10017 
International | 


i Se 


ene wpeoeeeeneeeeeeeaeoerereoerevenes 


